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1.0  INTRODUCTION 

To  enhance  the  understanding  of  flow  quality  in  the  AEDC  Supersonic  Wind  Tunnel  (A) 
and  Hypersonic  Wind  Tunnel  (B)  and  to  provide  documentation  of  the  wind  tunnel  test 
section  turbulence,  a  comprehensive  experimental  and  analytical  program  was  conducted 
during  the  period  October  1977  to  September  1980.  This  report  presents  a  comprehensive 
compilation  of  fluctuating  measurements  from  these  wind  tunnels. 

The  objectives  of  this  program  were  as  follows:  (1)  to  experimentally  evaluate  and 
document  the  test  section  turbulence  levels  in  Tunnels  A  and  B;  (2)  to  assess  any  differences 
found  between  dynamic  pressure  measurements  and  hot-wire  anemometer  results;  and  (3)  to 
formulate  and  implement  lest  procedures,  data  reduction  techniques,  and  software  needed 
for  obtaining  hot-wire  and  dynamic  pressure  data. 

The  scope  of  rhis  report  includes  the  identification  of  some  of  the  general  operational 
and  physical  features  of  the  facility  which  may  directly  or  indirectly  affect  the  tunnel 
turbulence  measurements  or  the  interpretations  of  dynamic  measurements  in  the  wind 
tunnel.  The  experimental  phase  covers  the  Mach  number  and  Reynolds  number  range  iti 
Tunnel  A  (M^  =  1 .5  to  5.0  and  RE  =  0.5  to  6  x  I06  per  foot)  and  in  Tunnel  B  (Mw  =  6.0 
and  8.0,  and  RE  =  0.5  to  5.0  x  10fi  per  foot).  The  measurements  include  dynamic  pressure 
data  obtained  from  surface-mounted  instrumentation  on  a  Hat  plate  and  a  sharp-nosed 
6-dcg  cone,  pitot  acousLic  probe  data,  and  hot-wire  anemometry  results.  The  report  also 
identifies  the  calibration  procedures  with  some  typical  calibration  results,  the  normalizing 
parameters  usually  applied  to  tunnel  wall  turbulence  measurements,  and  a  synopsis  of  the 
FFT  (Fast  Fourier  Transform)  signal  processing  used  in  this  study. 

2.0  BACKGROUND 

Wind  tunnel  test  section  turbulence,  its  intensity,  spatial  distribution,  and  spectral 
content,  has  become  increasingly  important  in  the  analysis  of  test  data.  The  sector- 
supported  model  in  ihe  lest  section  and  the  model  instrumentation,  particularly  the  dynamic 
pressure  sensors  (transducers)  and  the  hot-wire  survey  probes,  are  directly  affected  by  the 
test  section  turbulence.  Aerodynamically,  test  section  turbulence  can  affect  the  character  of- 
the  flow  over  the  model,  particularly  the  turbulence  or  instability  within  the  model  boundary 
layer  and  the  location  of  boundary-layer  transition  on  the  model.  Therefore,  it  is  important 
that  these  turbulence  and  noise  factors  be  measured,  identified,  and  documented  so  that  a 
proper  interpretation  and  correlation  can  be  made  of  flow  fluctuation  measurements 
oblained  in  the  wind  tunnel. 

In  a  high-speed  wind  tunnel  (M0CP>  >  1),  particularly  in  an  unheated  test  unit,  a  major 
source  of  fluctuating  disturbances  in  the  free-stream  flow  is  the  turbulent  tunnel  wall 
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boundary  layer  in  the  form  of  pressure  fluctuations  (a  sound  source),  as  noted  in  Refs.  I 
through  14.  Briefly,  this  tunnel  flow  disturbance  mechanism  is  viewed  as  a  moving  sound 
source  which  has  been  modeled  by  examining  the  characteristics  of  wavelets  produced  by  a 
turbulent  boundary  layer  in  a  supersonic  stream  (Refs.  1 1  and  12).  The  wavelets  produced  by 
the  shear  layer  arc  rationalized  by  assuming  that  a  moving,  fluid-dynamic  type  of  wavy  wall 
is  submerged  within  the  shear  layer.  This  wavy  wall  disturbance  moves  downstream  at  a 
velocity  greater  than  the  local  speed  of  sound.  The  wavy  wall  concept  supposedly  can  be 
attributed  to  the  velocity  or  vorticity  fluctuations  within  the  boundary  layer  which  in  turn 
are  viewed  as  internal  shear  layer  disturbances  produced  by  moving  “finite-shaped  bodies” 
or  “randomly  spaced  pulses”  (Ref.  1).  These  moving  fluid-dynamic  bodies  represent  finite¬ 
sized  cells  of  vorticity  or  circulation  moving  within  the  turbulent  shear  layer  which  transmit 
pressure  pulses  (wavelets)  into  the  exlernal  stream  producing  the  test  section  pressure 
fluctuations. 

This  concept  of  a  moving  sound  source  generated  by  the  turbulent  tunnel  wall  boundary 
layer  has  been  used  to  model  analytically  certain  tunnel  flow  fluctuating  properties.  With  a 
number  of  simplifying  assumptions  such  as  small  (relative  to  mean  flow  properties)  localized 
fluctuations  in  a  two-dimensional  flow  regime,  the  surface  pressure  fluctuation  is 
approximately  proportional  to  the  local  surface  shear  (Refs.  15,  16),  that  is,  P  is 
proportional  to  t*.  Therefore,  as  aptly  shown  in  Ref.  15,  the  rnis  (rooi-mean-square)  level 
of  the  pressure  fluctuations  can  be  described  in  terms  of  the  following  scmi-empirical 
relationship: 


~  f  2A  tO. 0296) /RE®'  2 1 

P'V=A.cr, - — - - - * -  (1) 

I0.5[l  +  (Tflr/T)l  -  0.22?;[(y  -  I)/2|M*I°-m 

Also,  the  power  spectral  distribution,  F(f).  of  the  stream  pressure  fluctuation  can  be 
shown  analytically  (Refs.  1  and  15)  and  by  well  defined  experimental  measurements  (Refs.  1, 
9,  10  and  1 1)  to  be  approximated  by  the  following  function: 

Fjj( 0/P2  =  F !  1 ,0/[  1  ^  (2ffl6*/uc)2| !  t2) 

This  simply  suggests  the  procedure  for  normalizing  tunnel  turbulence  spectral 
distributions;  it  further  suggests  that  the  spectral  distribution  approaches  asymptotically 
some  finite  upper  limit  as  the  frequency  (0  decreases  toward  the  d-c  level. 

Some  of  the  significant  wind  tunnel  and,  in  some  cases,  model  parameters  which  will  aid 
in  the  correlation  of  tunnel  turbulence  measurements  are  identified  and  discussed  in  the 
literature  (for  example  see  Refs.  1,  3,  8,  10,  and  15).  These  normalizing  parameters  include 
the  following: 
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1.  For  a  length,  the  boundary-layer  displacement  or  momentum  thickness. 

2.  For  time,  the  ratio  of  the  boundary-layer  thickness  parameter  to  the  free-stream 
or  local  sound  source  convective  velocity. 

3.  For  a  force  per  unit  area,  the  local  wall  shear  value  or  local  dynamic  pressure. 

4.  For  the  velocity,  the  convective  velocity  of  the  sound  source  or  the  free-stream 
velocity. 

These  are  the  major  parameters  that  will  be  considered  in  the  present  document. 

Previous  attempts  to  measure  flow-field  fluctuations  in  Tunnels  A  and  B  were  not  always 
successful.  As  noted  or  implied  in  Refs.  3,  5,  6,  and  7,  there  exists  a  fundamental 
disagreement  between  the  hot-wire  and  dynamic  pressure  (microphone)  results.  Donaldson 
(Ref.  7)  and  Laderman  (Ref.  3)  indicated  that  although  the  hot-wire  results  were  reasonable 
(that  is,  consistent  with  measurements  made  in  similar  facilities),  the  dynamic  pressure 
measurements  obtained  from  microphones  appear  to  be  unrealistically  high. 

Several  factors  are  known  to  affect  the  performance  or  response  of  a  dynamic  pressure- 
transducer.  In  the  case  of  a  microphone,  the  response  of  the  transducer  is  influenced  by  the 
natural  resonance  of  the  gage,  increasing  the  rms  of  the  gage  output.  Also,  the  scale  factor 
or  response  level  of  the  transducer  changes  with  the  local  temperature  and  pressure  level  of 
the  transducer  environment.  In  the  present  study,  an  effort  was  made  to  calibrate  the 
transducers  statically  and  dynamically  in  a  subatmospheric  environment  closely  approaching 
those  conditions  existing  in  the  Tunnel  A  and  B  environments.  In  general,  the  resonant 
frequency  of  most  of  the  transducers  used  in  this  study  were  well  above  the  maximum 
frequency  of  the  fluctuation  levels  sensed  in  the  tunnel  flow. 

Another  factor  which  affects  the  dynamic  pressure  measurement  is  associated  with  the 
transducer  installation  (Refs.  17  and  18).  Hanly’s  investigation  (Ref.  17)  conducted  on  a 
tunnel  wall  with  supersonic  flow  (Mach  numbers  1.68  and  2.5)  revealed  that  the  most 
consistent  and  reliable  results  were  obtained  when  the  transducer  was  flush  mounted  or 
slightly  recessed  (less  than  0.00055)  in  the  tunnel  (or  model)  wall.  Dynamic  pressure  gages 
protruding  as  little  as  0.05  percent  of  the  total  boundary-layer  thickness  (5)  produced 
significant  changes  in  the  gage  output  and  in  the  outputs  of  gages  located  immediately 
downstream  of  the  protruding  gage.  Increasing  the  local  stream  velocity  tended  to  increase 
the  gage  sensitivity  to  effects  produced  by  a  slightly  protruding  gage.  These  studies  (Refs.  17 
and  18)  demonstrate  that  the  transducer  protuberance  affects  the  mid-range  of  the  spectrum 
of  recorded  transducer  response.  In  the  present  study  every  effort  was  made  to  insure  (hat 
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wall  or  surface-mounted  transducers  were  either  flush  with  the  model  surface  or  slightly 
recessed  (about  0.005  in.). 

As  found  in  subsonic  tunnels  (Refs.  18  and  19)  and  in  the  supersonic  and  hypersonic 
wind  tunnels  (Refs.  5  and  7),  the  pressure  transducer  outputs  are  sensitive  to  acceleration 
effects,  electronic  noise,  and  sound  sources  transmitted  through  the  structural  members 
supporting  the  probe  model  containing  the  dynamic  pressure  sensors  (DPS).  One  way  to 
monitor  these  extraneous  effects  is  to  compare  the  output  of  an  active  gage  (one  exposed  to 
the  tunnel  flow)  with  the  output  from  a  dummy  or  inactive  gage  (one  isolated  from  the 
tunnel  flow  but  installed  immediately  adjacent  to  an  active  gage),  in  general,  these  external 
disturbances  tend  to  distort  (increase)  the  amplitude  of  the  pressure  signal  at  the  lower 
frequencies.  The  pressure  spectra  obtained  from  the  hot-wire  analyses,  the  more  reliable 
pressure  measurements  (see  Refs.  9,  10,  18,  19,  and  20),  and  theoretical  models  (see  for 
example  Refs.  1  and  15)  all  support  the  contention  that  the  pressure  fluctuation  spectrum 
should  decrease  abruptly  as  the  frequency  approaches  zero  (the  d-c  level).  However,  most  of 
the  pressure  spectra  obtained  from  wind  tunnel  measurements  where  the  lower  frequency 
components  of  the  spectra  were  retained  tend  to  increase  abruptly  at  the  lower  frequencies, 
thus  supporting  the  idea  that  “spurious”  pressure  fluctuations  (discussed  in  Ref.  19)  are 
being  sensed  by  the  transducer.  Therefore,  the  true  rms  level  of  a  pressure  signal  associated 
with  the  turbulence  in  the  flow  must  be  evaluated  with  these  spurious  components  of  the 
signal  spectrum  identified  and  removed  or  reduced  to  a  realistic  level. 

The  hot-wire  measurements  in  Tunnel  A  were  used  to  identify  pressure  fluctuations  in 
the  tunnel  test  section.  As  mentioned  previously  this  disturbance  is  introduced  primarily  by 
the  turbuleni  boundary  layer  along  the  tunnel  wall.  At  the  lower  Mach  numbers  (Mach 
number  <  2.0),  velocity  (vorticity)  fluctuations  may  be  introduced  by  the  stilling  chamber 
screens.  At  the  higher  Tree-stream  Mach  numbers,  these  screen  disturbances  are  usually 
damped  out  by  the  isentropic  compression  and  expansion  through  the  supersonic  nozzle. 
Temperature  (or  entropy)  fluctuations  are  generally  negligible  since  Tunnel  A  operates  at 
moderate  temperatures  (on  the  order  of  100  to  200° F).  In  Tunnel  B,  which  operates  at 
higher  stagnation  temperaiures,  particularly  at  Mach  number  8,  where  the  normal  operating 
stagnation  temperature  is  890°F,  the  temperature  fluctuations  or  spottiness  becomes 
measureable  and  must  be  taken  into  consideration  in  the  hot-wire  analysis.  As  indicated  in 
Ref.  2,  until  pressure  fluctuation  measurements  obtained  from  dynamic  pressure  sensors 
become  more  reliable,  hot-wire  anemometry  appears  to  be  the  most  reliable  means  for 
identifying  the  turbulence  in  a  supersonic/hypersonic  test  regime  (Tunnels  A,  B,  and  C). 

In  general,  the  measurements  made  in  Tunnel  A  were  obtained  with  50-pin. -diam  hot¬ 
wire  probes  and  in  Tunnel  B,  at  Mach  numbers  6  and  8,  with  20-pin. -diam  probes.  Although 
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the  acquisition  of  hot-wire  measurements  at  hypersonic  speeds  in  Tunnel  B  is  a  fairly 
standard  operational  procedure,  the  results  obtained  in  Tunnel  A  at  lower  supersonic  speeds 
are  unique  for  a  facility  of  this  size  {40-by  40-in.  test  section).  In  comparison  to  Tunnel  B, 
the  Tunnel  A  environment  is  much  more  severe  in  that  the  hot  wires  are  exposed  to  greater 
dynamic  loads,  greater  flow-field  particle  contamination,  and  higher  free-stream  humidity 
levels.  This  study  shows  that  in  spite  of  these  severe  environmental  effects  and  the  slightly 
higher  loss  of  hot  wires,  it  is  still  possible  to  obtain  useful  hot-wire  measurements  in  Tunnel 
A.  ■ 


3.0  TEST  FACILITIES 


3.1  TUNNEL  A 

3.1.1  Description 

Tunnel  A  (Fig.  1)  is  a  continuous,  closed-circuit,  variable  density  wind  tunnel  with  an 
automatically  driven  flexible-plate-type  nozzle  and  a  40-  by  40-in.  test  section.  The  tunnel 
can  be  operated  at  Mach  numbers  from  1.5  to  6  at  maximum  stagnation  pressures  from  29  to 
200  psia,  respectively,  and  at  stagnation  temperatures  up  io  750°R  at  Mach  number  6. 
Minimum  operating  pressures  range  from  about  1/10  to  1/20  of  the  maximum  at  each  Mach 
number.  The  tunnel  is  equipped  with  a  model  injection  system  which  allows  removal  of  the 
model  from  the  test  section  while  the  tunnel  remains  in  operation.  A  description  of  the 
tunnel  and  airflow  calibration  information  may  be  found  in  the  Test  Facilities  Handbook 
(Ref.  21). 

In  addition  to  the  turbulent  boundary  layer  on  the  nozzle  walls,  there  are  other  potential 
sources  of  disturbance  generated  in  the  Tunnel  A  circuit,  and  some  of  these  will  be  identified 
and  discussed.  Figure  lc  is  an  illustration  of  the  valves,  ducting,  screens,  and  stilling 
chamber  located  immediately  upstream  of  the  Tunnel  A  nozzle.  The  pressure  level  in  the 
tunnel  circuit  is  controlled  in  part  by  the  addition  of  high-pressure  air  through  valve  337  or 
237.  Each  valve  operates  by  opening  and  closing  at  regular  intervals  to  maintain  the  tunnel 
pressure.  At  low  Mach  numbers,  this  operation  generates  strong  pressure  pulses  which  can 
interfere  with  the  recording  of  dynamic  data.  Also,  the  addition  of  the  air  through  the  valve 
creates  a  sound  source  which  can  propagate  through  the  tunnel  structure  and  through  the 
flow  in  the  circuit.  For  testing  involving  dynamic  measurements,  a  manually  operated  valve 
(017)  located  downstream  of  the  diffuser  provides  a  more  favorable,  alternate  method  of 
adding  air  to  the  Tunnel  A  circuit. 

The  stilling  chamber  shown  in  Fig.  Id  contains  up  to  12  screens  of  various  mesh  sizes. 
The  purpose  of  the  screens  is  to  filter  the  air  flow  and  minimize  flow  gradients  entering  the 
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nozzle  region,  bul  at  the  lower  nozzle  Mach  numbers  these  screens  may  be  a  small  source  of 
velocity  fluctuations  which  may  reach  the  test  section  region  (see  Ref.  7). 

Downstream  of  the  tunnel  nozzle  region,  the  model  injection  tank  and  the  diffuser 
section  are  two  additional  sources  of  disturbance  which  may  affect  the  recording  of  dynamic 
measurements.  The  flow  over  the  slot  in  the  doors  over  the  model  injection  tank  is  a  source 
of  acoustical  disturbances  or  vibration.  The  turbulent  boundary-layer  interaction  with  the 
shock  wave  structure  in  the  diffuser  is  another  poreniial  source  of  disturbance. 

3.1.2  Test  Installation 

Dynamic  pressure  sensors  (DPS)  were  installed  at  various  stations  between  the  stilling 
chamber  and  test  section  of  Tunnel  A.  A  flush-mounted  DPS  unit  was  located  in  the  stilling 
chamber  as  shown  in  Figs,  lc,  d,  and  e.  As  noted  in  Fig.  le,  Lhc  0.2  in.-diam  transducer  was 
located  in  the  center  of  the  3.5-in.-diam  surface  with  nominally  a  0.25-in.  gap  between  the 
cylindrical  segment  (or  disk-shaped  surface)  and  the  stilling  chamber  wall.  Four  additional 
DPS  units  were  implanted  at  three  other  stations  in  the  wind  tunnel  nozzle;  two  DPS  units 
were  located  on  opposite  sides  of  the  tunnel  at  station  125  (see  Fig.  If),  one  at  station  30.7, 
and  one  in  the  model  injection  tank.  The  fixture  containing  the  DPS  units  inserted  into  the 
nozzle  side  walls  is  shown  in  Fig.  2.  A  pressure  tap  was  provided  so  that  the  local  nozzle  wall 
static  pressure  could  be  measured. 

The  first  test  installation  (1978)  for  the  study  consisted  of  a  mechanism  supporting  the 
flow-field  probes  which  could  be  traversed  relative  to  a  flat  plate  model  (Figs.  3  and  4).  The 
retractable  shielding  mechanism  shown  in  Fig.  3  for  the  hot-wire  probes  rotated  about  a 
fixed  point.  This  particular  design  was  subsequently  found  to  be  unsatisfactory  because  of 
the  strong  shock  system  produced  when  the  shield  was  retracted.  The  dynamic  pitot  probe 
was  mounted  on  top  of  the  retractable  probe-shielding  mechanism. 

In  the  1979  entry,  the  retractable  shielding  mechanism  was  modified  so  that  the  shield 
moved  almost  vertically  with  respect  to  the  probes  (see  Fig.  5).  Also,  a  new  pitot  acoustic 
probe  was  mounted  on  the  retractable  shield  fixture.  As  in  the  case  of  the  FY78  installation, 
the  probing  system  was  mounted  on  the  translating  mechanism  above  a  conical  body  which 
was  also  instrumented  with  DPS  units  as  shown  in  Fig.  6.  Both  the  FY78  and  FY79 
installations  could  be  retracted  into  the  model  injection  tank  for  additional  model  work  and 
replacement  of  hot-wire  probes. 
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3.2  TUNNNKL  B 

3.2.1  Description 

Tunnel  B  (Fig.  7)  is  a  closed-circuit  hypersonic  wind  tunnel  with  a  50-in. -diam  test 
section.  Two  axisymmetric  contoured  nozzles  are  available  to  provide  Mach  numbers  of  6 
and  8,  and  the  tunnel  may  be  operated  continuously  at  pressure  level  ranges  from  20  to  300 
psia  at  Mach  number  6,  and  from  50  to  900  psia  at  Mach  number  8,  with  air  supplied  by  the 
VKF  main  compressor  plant.  Stagnation  temperatures  sufficient  to  avoid  air  liquefaction  in 
the  test  section  (up  to850°Rat  Mach  No.  6  and  up  to  l,350°Rat  Mach  No.  8)  are  obtained 
through  the  use  of  a  natural  gas-fired  combustion  heater  located  250  ft  from  the  Tunnel  B 
nozzle.  The  entire  tunnel  (throat,  nozzle,  test  section,  and  diffuser)  is  cooled  by  integral, 
external  water  jackets.  The  tunnel  is  equipped  with  a  model  injection  system,  which  allows 
removal  of  the  model  from  the  test  section  while  the  tunnel  remains  in  operation,  A 
description  of  the  tunnel  may  be  found  in  the  Test  Facilities  Handbook  (Ref.  21). 

Except  for  the  Thermal  disturbances,  the  disturbances  produced  in  the  stilling  chamber  of 
Tunnel  B  are  attenuated  by  the  isentropic  compression  through  the  nozzle.  As  in  the  case  of 
Tunnel  A,  (he  flow  over  the  Tunnel  B  model  injection  tank  and  the  unstable  How  in  the 
diffuser  section  aft  of  the  tunnel  test  section  are  additional  sources  for  disturbances  that  may 
be  sensed  by  dynamic  pressure  probes. 

3.2.2  Test  Installation 

For  the  1978  testing  the  flat  plate  used  in  the  Tunnel  A  studies  was  mounted  on  the 
Tunnel  B  sector.  The  probes  were  supported  on  the  x-z  probe  mechanism  shown 
schematically  in  Fig.  8. 

In  1979  the  tests  used  the  x-z  probing  mechanism  and  the  sector-mounted  6-deg  sharp 
cone  as  shown  in  Fig.  8a.  A  pitot  acoustic  probe  was  strut  mounted  to  the  sting  of  the  sharp 
cone  model.  The  pitot-acoustic  probe  was  removed  From  the  installation  while  flow-field 
surveys  were  made  with  [he  x-z  probe  mechanisms.  The  pitot  acoustic  probe  could  not  be 
exposed  to  the  high  thermal  environment  of  Tunnel  B  for  the  longer  periods  of  time  required 
lo  obtain  flow-field  surveys  with  (he  x-z  probe  mechanism. 
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3.3  PLAINT 


The  VKF  main  compressor  plant  operates  over  an  extremely  wide  range  of  pressures, 
pressure  ratios,  and  mass  flow  rates  to  meet  the  requirements  for  closed  circuit  operation  of 
the  continuous  How  Tunnels  A,  B,  and  C.  A  detailed  description  of  this  plant  is  given  in  Ref. 
21 .  The  plant  machinery  (the  drive  motors  and  compressors  —  see  Figs.  2. 1  and  2.2  of  Ref. 
21)  generates  vibrations  close  to  the  rotational  speeds  and  harmonics  of  these  speeds  of  the 
various  stages  of  plant  operation.  These  vibrations  are  transmitted  through  the  structures  to 
the  Tunnel  A,  B,  and  C  tunnel  areas.  A  summary  of  the  rotational  speeds  is  presented 
below. 


Stages  of 
Operation 

Plant  Operation 
Nominal 

Rotational  Speeds  (Hz) 

Tunnel  Mach  Number 

1 

60 

1.5  to  2.0  (Tunnel  A.) 

2 

60 

1.5  to  3.0 

3 

47 

3.0  to  4.0 

4 

68 

4.0  to  5.5 

5 

99 

5.5  to  6.0 

6 

162 

5.5  to  6.0 

7 

201 

6.0,  8.0  (Tunnel  B) 

6 

203 

10.0  (Open  Circuit) 

9 

200 

10.0 

For  example,  on  the  shaft  of  compressor  81  (which  is  included  in  the  eighth  stage)  a 
significant  vibration  at  60  Hz  with  secondary  frequencies  of  20,  40,  124,  and  220  Hz  was 
recorded.  Significant  vibrations  at  108  and  208  Hz  were  recorded  on  the  shaft  of  compressor 
71  when  it  was  operating  in  an  unusual  manner.  Of  course,  several  stages  of  the  plant  were 
operating  when  the  aforementioned  vibrations  were  detected  and  were  simply  transmitted 
to,  for  example,  compressors  81  and  71.  These  vibrations  have  been  identified  to  explain 
potential  acceleration  modes  which  might  be  transmitted  and  detected  or  superimposed  on 
dynamic  signals  recorded  in  the  tunnel  test  section  area. 


4.0  APPARATUS 


4.1  MODELS 

A  flat  plate  model  (in  1978)  and  a  sharp  cone  model  (in  1979)  were  used  to  measure  wall 
static  pressure  fluctuations  on  representative  model  surfaces  in  the  test  sections  of  Tunnels  A 
and  B.  These  model  measurements  were  used  in  a  comparison  with  the  static  pressure 
fluctuations  inferred  from  the  hot-wire  anemometer  measurements  and  the  dynamic  pitot 
pressure  measurements  made  during  the  present  investigation. 
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4.1.1  Flal  l’lale  Model 

The  flat  plate  model  (Fig.  9)  had  a  length  oT  10  in.  and  a  span  of  8  in.  and  was  fabricated 
of  stainless  steel.  The  plate  was  instrumented  with  two  flush-mounted  piezoelectric  dynamic 
pressure  transducers,  iwo  static  pressure  taps,  two  thermocouples,  and  three  flush-mounted 
thermopile  heat-transfer  gages.  The  dynamic  pressure  transducers  had  circular  sensing  areas 
with  nominal  diameters  of  0.2  and  0.4  in.  that  were  centered  at  stations  5  in.  from  the 
leading  edge  of  the  flat  plate  and  0.5  in.  on  either  side  of  the  axial  centerline  of  the  plate.  The 
transducers  were  mounted  in  nylon  inserts  in  the  flat  surface  of  the  model  and  were  installed 
as  indicated  in  Fig.  9.  The  use  of  this  method  of  mounting  the  transducers  was  subsequently 
judged  to  be  unsatisfactory,  as  discussed  in  Section  4.2.2.  The  pressure  orifices,  which  had 
an  inside  diameter  of  0.064  in.,  were  spaced  2  in.  apart,  in  the  axial  direction.  The  static 
pressures  were  used  to  confirm  alignment  of  the  flat  surface  with  the  free-stream  flow.  The 
thermocouples  were  Chromel®  -Alumel®  junctions  finished  flush  with  the  flat  surface  and 
were  used  primarily  to  monitor  the  temperature  of  the  model  to  safeguard  the  dynamic 
pressure  sensors.  The  heat-transfer  gages  were  included  with  the  instrumentation  in 
anticipation  of  the  use  of  the  model  with  roughness  elements  applied  to  the  surface  to 
promote  boundary-layer  transition  over  the  dynamic  pressure  sensors.  However,  time  did 
not  permit  this  extension  of  the  study.  The  model  was  equipped  with  sideplates  to  contain 
the  higher-pressure  air  on  the  underside  of  the  model.  A  shadowgram  of  the  flat  plate  model 
in  the  Mach  number  8  flow  in  Tunnel  B  is  presented  in  Fig.  10. 

4.1.2  Sharp  Cone  Model 

The  sharp  cone  model  (Fig.  1 1)  was  a  6-deg  (half-angle)  right-ahgle  cone  with  a  length  of 
approximately  43  in.  and  a  base  diameter  of  9  in.  and  was  fabricated  of  stainless  steel.  The 
cone  was  instrumented  with  dynamic  pressure  transducers,  static  pressure  taps,  and 
thermocouples.  Surface  pressure  fluctuations  were  measured  using  nine  dynamic  pressure 
transducers  located  at  three  axial  stations  along  the  model.  The  sensor  mourning  stations 
were  chosen  so  that  the  transducers  at  the  forward  and  aft  positions  would  be  located  under 
laminar  and  turbulent  boundary-layer  nows,  respectively.  The  cone  boundary  layer  over  the 
transducers  mounted  at  the  central  station  was  laminar,  transitional,  or  turbulent  depending 
upon  tunnel  test  conditions.  In  addition  one  dynamic  pressure  transducer  was  used  as  a 
“blind”  transducer  in  the  cone  to  assess  the  acceleration  effects  imposed  upon  the  surface- 
mounted  sensors  by  model  vibration.  The  dynamic  pressure  transducers  installed  in  the  cone 
surface  were  mounted  in  molded  silicone  rubber  inserts  of  RTV®  -60  (red)  as  illustrated  in 
Fig.  1  lb.  The  insert  partially  isolated  the  sensor  from  model  vibrations  and  thereby  reduced 
associated  acceleration  effects.  For  the  measurements  in  Tunnel  B  a  0.005-in.-thick  cap  of 
RTV-732  (white)  was  added  to  cover  the  transducer  sensing  surface  and  provide  protection 
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from  the  high-temperature  environment.  (The  white  compound  has  better  adhesive 
characteristics  than  the  red  compound,  according  to  manufacturer’s  specifications.)  Each 
cap  was  contoured  with  the  cone  surface;  however,  the  small  diameters  of  the  sensing  areas 
of  the  PCB®  -105  transducers  resulted  in  only  a  small  indentation  in  the  model  contour 
when  caps  were  not  used  in  Tunnel  A.  The  blind  transducer  was  threaded  into  a  metal  insert 
implanted  in  the  inner  surface  of  the  model  wall  (Fig.  11)  and  operated  in  the  environment, 
particularly  the  pressure,  of  the  interior  cavity  of  the  cone  model. 

Eight  static  pressure  measurements  were  obtained  using  surface  pressure  taps  located 
circumferentially  at  90-deg  intervals  around  the  model  at  two  stations  near  the  forward  and 
central  transducer  mounting  stations.  The  static  pressure  measurements  provided  a  direct 
indication  of  the  cone  alignment  with  the  free-siream  flow  direction  and  were  used  in 
reducing  boundary-layer  survey  data. 

In  addition  to  the  rubber  caps  discussed  earlier,  the  following  provisions  were  made  to 
protect  the  dynamic  pressure  transducers  in  the  high-temperature  environment  ot  Tunnel  B: 
the  nose  of  the  model  was  integrally  water-cooled  when  the  model  was  exposed  to  the  tunnel 
flow  for  long  periods  of  time,  and  six  thermocouples  were  used  to  monitor  model  wall 
temperature.  Three  Chromel-AIumel  thermocouples  were  attached  to  the  inside  wall  of  the 
cone,  and  three  Chromel-eonstantan  coaxial  thermocouples  were  mounted  flush  with  the 
exterior  surface  of  the  model.  Sensor  locations  are  indicated  in  Fig.  1 1. 

4.2  Probes 

The  principal  measurements  in  the  frec-stream  flow1  in  the  test  section  ol  each  tunnel 
were  made  using  hot-wire  anemometer  and  dynamic  pitot-pressure  probes.  These 
measurements  were  supplemented  using  conventional  pitot  pressure  and  total  temperature 
probes.  All  probes  were  fabricated  by  personnel  at  AEDC. 

4.2.1  Anemometer  Probes 

The  hot-wire  anemometer  probes  were  made  using  platinum  -  10  percent  rhodium  wires, 
drawn  by  the  Wollaston  process,  of  20-' and  50-/xin.  nominal  diameter  and  approximately 
150  diameters  length.  The  sensor  wire  was  attached  to  sharpened  3-mil  nickel  wire  supports 
using  a  gold-bonding  technique  developed  by  Philco-Ford  Corporation  (Ref.  22).  The  probe 
body  (Fig.  12)  was  composed  of  an  alumina  cylinder  of  0.031-in.  diam  and  0.40-in.  length 
bonded  to  a  second  alumina  cylinder  of  0.093-in.  diam  and  3.0-in.  length  using  ceramic 
cement.  The  cylinders  (thermocouple  insulators)  were  manufactured  with  two  holes  along 
their  length  that  were  used  for  the  leads  of  the  sensor.  The  nickel  wire  supports  that  passed 
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through  the  smaller  cylinder  were  silver-soldered  to  copper  leads  that  were  inserted  into  the 
larger  cylinder.  Ceramic  cement  was  used  at  each  end  of  the  probe  body  to  secure  the 
supports  and  leads.  The  physical  and  heat-transfer  characteristics  of  each  hot-wire 
anemometer  probe  used  for  measurements  were  determined  by  calibration  procedures 
discussed  in  Appendix  A.  The  probes,  which  were  designed  to  meet  the  structural  and  high- 
frequency  requirements  for  measurements  in  hypersonic  flow,  performed  well  at  Mach 
.numbers  as  low  as  3.  Only  a  very  limited  amount  of  anemometer  data  was  acquired  at  Mach 
numbers  1 .5  and  2.  It  is  believed  that  the  principal  causes  of  wire  loss  at  the  lower  Mach 
numbers  were  particulate  matter  in  the  flow  and  the  high  aerodynamic  loads  on  the  wire. 

4.2.2  Dynamic  Pitot  Pressure  Probes 

Dynamic  pitot  pressure  probes  of  two  designs  (Fig.  13)  were  used  in  this  investigation.  In 
the  initial  design  a  single  dynamic  pressure  transducer  was  mounted  in  a  frustum-shaped 
stainless  steel  tip  threaded  to  hold  the  transducer  with  no  provision  for  isolation  from  the 
probe  body  vibration  (Fig.  13a).  The  mounting  included  a  thin  (of  the  order  of  0.015  in. 
thick)  silicone  rubber  (RTV-732)  cap  to  protect  the  transducer  sensing  surface  from  particles 
in  the  flow  in  the  tunnels  and  the  more  intense  Tunnel  B  thermal  environment. 

The  second  probe  design  (Fig.  13b)  was  developed  to  furnish  some  degree  of  isolation  of 
the  transducer  from  probe  body  vibration  effects  and  included  provision  for  mounting  an 
identical  transducer  to  serve  as  a  blind  sensor  to  assess  the  remaining  vibration  effects.  The 
two  transducers  were  mounted  in  molded  silicone  rubber  inserts  which  were  separated  by  a 
nylon  sleeve.  The  assembly  was  held  in  the  probe  body  by  a  threaded  plug  at  the  base  of  the 
probe  and  was  separated  from  rhe  rubber  insert  holding  the  blind  transducer  by  a  second 
nylon  sleeve.  The  measurements  made  using  the  second  probe  design  were  found  to  be 
characterized  by  a  higher  signal-to-noise  ratio  than  those  measurements  made  with  the  probe 
of  the  initial  design.  The  sensing  surface  of  the  primary  transducer  was  protected  by  a  thin 
silicone  rubber  (RTV-732)  cap,  as  in  the  first  probe  design.  It  should  be  noted  that  the 
dynamic  pressure  transducers  used  in  these  probes  may  be  unstable  when  the  temperature  of 
the  transducer  exceeds  250°F.  As  a  consequence,  the  length  of  time  that  the  probe  was 
continuously  exposed  to  the  free-stream  flow  in  Tunnel  B  was  quite  short  (less  than  one 
minute).  f 

4.2.3  Pitot  Pressure  and  Total  Temperature  Probes 

The  pitot  pressure  probe  (Fig.  14)  was  designed  for  boundary-layer  surveys  on  the 
models  in  the  vicinity  of  selected  surface-mounted  dynamic  pressure  transducers.  The  probe 
was  made  by  flaiiening  a  0.032-in.  O.D.  (0. 022-in.  I.D.)  tube  to  an  inside  vertical  height  of 
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0.005  in.  for  the  probe  orifice.  The  tube  was  bent  to  align  it  parallel  to  the  model  surface 
over  a  length  of  0.38  in.  from  the  probe  orifice.  The  probe  was  also  used  to  determine  local 
Mach  number  in  the  vicinity  of  the  anemometer  probe  tip. 

The  total  temperature  probe  (Fig.  14)  was  fashioned  from  a  length  of  sheathed 
thermocouple  leads  (Chromel-Alumel)  with  wires  of  1 .5-mil  diam.  The  wires  were  bared  for 
a  length  of  15  mils,  and  a  thermocouple  junction  of  5-mil  diam  was  formed  at  the  end  of  the 
wires.  The  thermocouple  was  used  without  shielding  to  reduce  the  probe  dimensions  and 
make  it  more  comparable  to  those  of  the  pitot  pressure  probe.  The  probe  was  used  in  Tunnel 
B  for  free-stream  and  boundary-layer  measurements. 

4.2.4  Probe  Holders 

Holders  for  mounting  the  anemometer,  pitot  pressure,  and  total  temperature  probes  as 
required  were  fabricated  with  0.1 25-in. -diam  (O.D.)  tubes  in  which  the  position  of  each 
probe  could  be  individually  adjusted  so  that  the  sensing  tips  or  orifice  would  be  located  at 
the  same  model  axial  station  for  surveys  of  the  model  boundary  layer.  Centerlines  of  the 
tubes  of  the  probe  holders  were  laterally  spaced  approximately  0.13  in.  apart  in  a  common 
plane.  Each  probe  was  held  in  place  by  a  set  screw,  facilitating  efficient  replacement  of 
probes.  The  holder  used  in  Tunnel  A  measurements  accommodated  two  probes,  and  the 
holder  of  the  ihree  probes  used  in  Tunnel  B  could  accommodate  as  many  as  four  probes. 

4.3  INSTRUMENTATION 

4.3.1  Hot-Wire  Anemometer 

Constant-current  hot-wire  anemometer  instrumentation  was  used  for  the  present 
investigation.  The  anemometer  current  control  (Philco-Ford  Model  ADP-13),  which 
supplies  the  heating  current  to  the  sensor  wire,  is  capable  of  maintaining  the  current  at  any 
one  of  15  preset  (adjustable)  levels  individually  selected  using  pushbutton  switches.  The 
anemometer  amplifier  (Philco-Ford  Model  ADP-12),  which  amplifies  the  wire  response 
(a-c  voltage)  signal,  contains  the  circuits  required  to  electronically  compensate  the  signal  for 
the  thermal  lag  that  is  a  result  of  the  finite  heat  capacity  of  the  wire.  A  square-wave 
generator  (Shapiro/Edwards  Model  G-50)  was  used  to  determine  the  time  constant  of  the 
sensor  wire  whenever  required.  The  wire  heating  current  and  mean  voltage  were  fed  to 
autoranging  digital  voltmeters  for  a  visual  display  and  to  a  magnetic  tape  machine  for 
recording  (see  Section  4.3.3).  The  sensor  response  (a-c  voltage)  was  fed  to  an  oscilloscope 
for  visual  display  of  the  real-time  signal  and  to  a  wave  analyzer  (Hewlett-Packard  Model 
8553B/8552B)  for  visual  display  of  the  spectrum  of  the  fluctuating  signal  in  real  time;  the 
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sensor  response  was  recorded  on  magnelic  tape  for  subsequent  analysis.  A  detailed 
description  of  the  VKF  hot-wire  anemometer  instrumentation  is  given  in  Ref.  23.  In  addition 
to  the  analog  data  records,  a  digital  record  of  the  wire  heating  current,  mean  voltage,  and 
rms  voltage  of  Lhe  sensor  response  is  acquired  (as  of  1979)  through  the  appropriate  VKF 
tunnel  data  system. 

4.3.2  Dynamic  Pressure  Sensors 

Two  types  of  transducers  were  used  to  record  the  sound  spectra  existing  in  the  wind 
tunnel.  The  primary  measuring  instruments  were  piezoelectric  sensors  with  diameters  of 
0.099-,  0.21 8-,  and  0.405  in.  These  are,  basically,  quartz  transducers  in  which  the  quartz 
element  produces  a  change  in  response  to  an  abrupt  change  in  the  applied  pressure.  Special 
calibration  procedures  were  used  in  determining  the  dynamic  response  of  these  transducers 
(Appendix  B).  Depending  on  the  size  of  the  transducer,  the  sensitivity  generally  varied  from 
25  to  100  mv  per  psi,  the  resonant  frequency  varied  from  250  to  500  kHz,  and  the  rise  time  of 
the  transducers  was  nominally  1  to  2  /isec. 

The  remaining  transducers  used  in  this  experiment  were  resistance-type  sensors  which 
had  sensitivity  factors  of  nominally  12  to  16  mv  per  psi  and  a  frequency  range  from  1/5  to 
1/10  of  the  range  of  the  piezoelect ronic  sensors.  The  uncertainty  in  the  dynamic  pressure 
measurements  is  based  on  significantly  fewer  calibration  points  (degrees  of  freedom); 
therefore,  a  higher  Student’s  “t”  value  must  be  applied  to  the  calibration  standard 
deviation,  and  this  indicates  a  higher  measurement  uncertainty.  This  is  not  an  “end-to-end” 
calibration  in  that  the  input  to  the  installed  sensing  gage  could  not  be  compared  directly  to 
the  final  displayed  output.  The  same  recording  and  signal  processing  instrumentation  used 
in  the  calibrations  was  used  in  the  wind  tunnel  tests. 

4.3.3  Signal-Processing  Equipment 

Conditioned  analog  response  signals  from  the  hot-wire  anemometer  and  dynamic 
pressure  transducers  were  recorded  on  the  VKF  Bell  and  Howell  model  VR3700B  magnetic 
tape  machine  in  ihe  FM  mode.  Each  channel  was  calibrated  and  adjusted  to  have  a  signal-to- 
noise  ratio  of  35  db  for  a  1.00-v  rms  output.  The  tape  machine  frequency  response  was  +  1 
to  -  3  db  over  a  frequency  range  from  dc  to  500  kHz.  In  the  present  calibration,  a  sine  wave 
generator  was  used  to  check  each  channel  at  several  discrete  frequencies,  using  an  rms 
voltmeter  which  was  periodically  checked  on  1-,  10-,  and  100-v  ranges  and  found  to  have  an 
uncertainty  of  less  than  ±  1.0  percent  of  the  reading  within  the  frequency  band  from  25  Hz 
to  100  kHz.  All  magnetic  tape  recordings  were  made  at  a  tape  speed  of  120  in. /sec. 

,  ; 
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Sensor  response  signals  recorded  on  magnetic  tape  during  hot-wire  anemometer  and 
dynamic  pressure  transducer  measurements  were  subsequently  processed  using  the  VKF 
spectrum  analyzer,  a  Spectral  Dynamics  Corporation  Model  SD360  digital  signal  processor. 
This  processor  contains  a  dedicated  computer  programmed  with  a  Fast  Fourier  Transform 
(FFT)  algorithm  to  obtain  single  channel  1,024-point  (or  dual  channel  5 12-point /channel) 
forward  or  inverse  Fourier  transforms.  Other  operations  often  required  in  signal  analysis 
(for  example,  calculations  of  auto-correlation  and  cross-correlation  functions  and  their 
respective  Fourier  transforms,  the  spectral  density  and  cross-spcciral  density)  can  be 
performed  by  the  processor.  The  frequency  range  for  an  analysis  can  be  chosen  from  among 
57  available  ranges,  full  scale  from  10  Hz  to  150  kHz,  and  each  range  has  a  pair  of  anti¬ 
aliasing  filters  with  120  db/octave  rolloff.  The  processor  uses  block  floating-point 
arithmetic  in  the  FFT  computations  to  avoid  loss  of  data  by  overflow.  Rectangular  and 
Kaiser-Bessel  weighting  functions  are  available  for  use  in  the  data  processing.  Linear 
ensemble  averaging  and  exponential  averaging  are  available  to  obtain  signal-to-noise 
enhancement  or  to  improve  the  statistical  confidence  in  the  results. 

The  SD360  analyzer  is  interfaced  with  a  VKF  Digital  Equipment  Corporation  PDP  1 1/40 
minicomputer  for  disk  storage  of  analysis  results  and  subsequent  transmission  of  the  results 
to  the  VKF  Digital  Equipment  Corporation  DEC  System  10  computer  for  further 
processing. 

5;0  DATA  ANALYSIS 

5.1  SIGNAL  PROCESSING  AND  NORMALIZATION 

5.1.1  Procedure 

As  noted  previously  (Section  4.3),  the  hot-wire  anemometry  and  dynamic  pressure 
signals  were  initially  recorded  on  magnetic  tape  for  further  analysis.  The  Spectral  Dynamics 
SD360  digital  signal  processor  (see  Section  4.3.3)  was  used  to  convert  the  recorded  real  time 
analog  signal  into  a  digitized  spectrum  in  the  frequency  domain  for  further  analysis  and 
normalization.  The  SD360  samples  an  input  signal  at  a  rate  equal  to  the  reciprocal  of  the 
product  of  2.048  times  the  selected  bandwidth  of  the  analysis. 

Since  the  playback  speed  of  the  analog  tape  can  be  varied  relative  to  the  tape  speed  of  the 
signal  recording,  the  true  bandwidth  of  the  signal  analysis  using  the  SD360  can  be  changed 
significantly.  For  example,  in  the  present  analysis,  the  Tunnel  B  hot-wire  results  were 
effectively  analyzed  over  a  total  signal  bandwidth  (spectrum)  of  600  kHz.  (The  playback 
speed  of  the  tape  machine  was  one-fourth  the  recording  speed,  and  the  analyzer  bandwidth 
was  set  at  150  kHz).  In  Tunnel  A  the  total  bandwidth  of  the  hot  wire  results  was  160  kHz 
because  the  tape  speed  ratio  was  one-half  and  the  analyzer  bandwidth  was  80  kHz. 
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The  DPS  results  were  processed  at  selected  total  signal  bandwidths  of  50  kHz,  100  kHz, 
and  150  kHz.  In  this  instance  the  playback  and  recording  speeds  were  identical;  that  is,  the 
tape  speed  ratio  was  one. 

The  SD360  was  also  used  to  digitize  real  time  analog  signals  such  as  the  real  time 
response  of  the  DPS  units  located  in  the  transient  flow  of  a  shock-tube  calibrator.  In  these 
situations,  the  bandwidth  of  the  sampling  rate  was  set  at  150  kHz,  and  analog  tape  playback 
speed  was  one-fourth  the  recording  speed  (see  Appendix  B  for  further  discussion).  The 
actual  real  time  signal  was  digitized  at  a  rate  of  0.8138  /isec  per  point.  The  total  length  of  the 
signal  sample  was  0.833  msec. 

5.1.2  Statistical  Uncertainties 


The  normalized  statistically  defined  error  (uncertainty)  associated  with  the  narrow-band 
spectral  analysis  follows  (Ref.  24). 

The  random  error  in  the  mean  square  value  of  the  signal  is  e(fms2)/fms2  *  (BW  ■  T)-°-5, 
or  simply  (N/2)-°  since  the  sample  length,  T,  equals  N/(2.048  BW).  The  term  "N”  is  the 
total  number  of  points  in  the  interval  T  used  to  define  rms2,  the  mean  squared  value  of  the 
signal.  For  the  SD360,  the  value  of  N  is  either  1,024  or  2,048;  therefore,  the  corresponding 
random  error  is  either  ±4.4  percent  or  ±3.1  percent. 

The  random  error  associated  with  the  amplitude  of  the  power  spectral  density  (PSD) 
distribution  is  defined,  according  to  Ref.  24,  as  e(F)/F  *  [Af(nd  •  T)]-°4.  since  Af  =  1/T, 
then  e(F)/F  =  nd-0-5.  The  SD360  can  take  up  to  21S  signal  samples  (nd)  and  average  the 
results.  A  summary  of  the  random  error  analyses  associated  with  the  use  of  the  SD360 
analyzer  is  given  in  the  following  table. 

ESTIMATED  RANDOM  ERRORS 


Item 

BW, 

kHz 

N 

e  (rma^) /rras^, 
percent 

nd 

e(F)/P, 

percent 

Hot  Wire  (Tunnel  A) 

160 

1,024 

±4.4 

26 

±12.5 

Hot  Wire  (Tunnel  B) 

600 

1 

1 

DPS  (Tunnel  A) 

100 

1 

DPS  (Tunnel  B) 

150 

? 

DPS  (Tunnel  A/B) 

50 

25 

±17.7 

5.1.3  Data  Processing 


The  Fourier  transform  provides  a  convenient  vehicle  for  transforming  real  lime  signals 
into  the  frequency  domain.  These  hol-wire  and  DPS  signals  are  analyzed  as  if  they  were  the 
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result  of  an  essentially  slationary  (ergodic)  random  process.  The  continuous  form  of  the 
transform  (Ref.  25)  is 

ec 

H<f)=/  h(t)e-J27rfl  dt  ,  units -sec  <3> 

•~CC 

and  the  inverse  transform  is 

oa 

h(l)=J  ll(ne-,2,rfldf  units  (4) 

— Os> 

Since  the  SD360  is  an  analog-to-digital  analyzer  with  an  FFT  processor,  the  character  of 
the  discrete  transform  takes  the  following  form: 

Hn  (^-)=--nf  h(mc-i2^/N)  ,  units  (5) 

The  scaling  between  the  discrete  transform  and  the  corresponding  amplitude  of  the 
transform  is  as  follows  (see  Ref.  26,  p.  102): 

H  =  *  rHHD f units  —  sec  (6) 

Where  “NT,”  the  total  length  of  time  of  the  sampled  signal,  equals  N/(2.048  •  BW)  and  the 
incremental  bandwidth  of  the  digital  spectral  analysis, 

AT  =  CN  -  T)~]  ,  sec-’  (7) 

The  power  spectral  density  is  normally  defined  as 

F(f)  =  lim  '  ■  f  li2(t;Af.f)  dt  .  units2  —  sec  ^ 

AlUo  ir- 1  Ja 

whereas  the  digital  form  is 

F|,  (  =  FD  (n  •  Af)  =  ~  [HD(n,A0  •  H*(r.,A0]  ,  units2/Sec  (9> 

where  Ho(n,  AO  is  the  conjugate  of  H i>(n,  AO-  The  corresponding  amplitudes  of  PSD  for 
the  continuous  function  are 


F  (n  ■  Af)  =  NT  |H  ()(n  •  Af)  1 2  ,  units2 -sec 


(10) 
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Since  the  root-mean-square  (rms2)  of  the  signal  is  normally 


2 

rms 


Rf)df 


units2 


(H) 


then  the  digital  form  becomes 


n/2 


tJ-L  |H(n-Af)|2 


n-1 


units2 


(12) 


This  summarizes  the  definitions  of  parameters  used  to  process  the  hot-wire  and  dynamic 
pressure  signals  where  the  “unit,”  the  amplitude  of  the  signal,  is  usually  in  millivolts  (mv). 

5.1.4  Normalization 

The  stilling  chamber  and  test  section  properties  and  probe  measurements  included  in  this 
study  are  listed  in  Appendix  C,  Table  1.  This  listing  contains  most  of  the  test  conditions 
covered  in  this  flow  fluctuation  study.  Some  of  the  tunnel  circuit  parameters  such  as  PT,  P, 
Pp,  and  q  were  used  in  the  normalization  of  the  fluctuating  parameters  and  also  in  the  other 
aspects  of  the  hot-wire  analyses.  A  description  of  the  normalization  procedure  and 
determination  of  normalization  parameters  is  given  in  Appendixes  C  and  D. 

5.2  ANEMOMETER  SIGNAL  ANALYSIS 

Hot-wire  anemometry  techniques  were  used  to  determine  the  levels  of  fluctuations  in  the 
free-stream  flow  along  the  centerline  of  the  test  sections  of  Tunnels  A  and  B.  The 
measurements  included  the  hot-wire  heating  current,  the  d-c  voltage  drop  across  the  wire, 
and  the  wire  a-c  voltage  response  to  fluctuations  in  the  flow.  The  a-c  voltage  component  was 
recorded  in  analog  form  on  magnetic  tape  for  subsequent  spectral  analysis,  and  the  rms  level 
of  the  a-c  voltage,  unfiltered  and  including  the  electronic  noise,  was  determined 
simultaneously  using  a  true-rms  meter.  The  procedure  used  to  process  the  data  and  to 
interpret  the  results  are  discussed  below. 

5.2.1  Spectral  Analysis 

The  anemometer  a-c  voltage  response  signals  recorded  on  magnetic  tape  were  used  as 
input  to  the  VKF  SD-360  spectrum  analyzer  (Section  4.3)  to  obtain  the  digitized  power 
density  spectrum  of  the  signals.  The  digital  output  of  the  SD-360  processor  consisted  of 
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1,024  points  equally  spaced  over  ihe  selected  frequency  range  of  the  analyzer.  (See  the  listing 
in  Section  5.1.2).  Approximately  830  spectra  were  obtained  from  the  free-stream 
measurements  made  in  the  investigation.  Each  spectrum  corresponded  to  one  wire  heating 
current  at  one  unit  Reynolds  number  at  one  Mach  number.  Repeat  measurements  were  made 
in  many  cases  generally  because  of  wire  failure  during  a  series  of  measurements.  Although 
all  of  the  spectra  were  examined,  only  approximately  30  percent  of  the  data  were  used.  In 
general,  data  were  rejected  for  the  present  investigation  under  the  following  conditions: 
(1)  if  they  were  part  of  an  incomplete  set  of  measurements;  (2)  if  the  power  density  spectrum 
of  the  wire  response  signal  was  not  clearly  above  that  of  the  electronic  noise  over  a 
“reasonable”  bandwidth;  (3)  if  spurious  signals  such  as  those  caused  by  wire  strain-gage 
effects  were  present  in  the  spectrum;  (4)  if  the  general  “shape”  of  the  specrrum  was 
inconsistent  with  that  of  the  other  data  of  the  same  set  of  measurements;  and  (5)  if  the  hot¬ 
wire  probe  experienced  a  shift  in  unheated-wire  resistance  which  was  unrelated  to  a  change 
in  tunnel  operating  conditions. 

5.2.2  Signal  Enhancement 

In  the  acquisition  of  data,  the  electronic  noise  of  the  assembled  instruments  of  the 
anemometer  system  was  recorded  at  each  test  condition  for  which  flow  fluctuation 
measurements  were  made.  The  noise  (unheated  wire)  recordings  were  processed  in  the  same 
manner  as  the  response  signals  of  the  heated  wire.  The  value  of  the  power  density  spectrum 
of  the  noise  was  subtracted,  at  each  discrete  frequency,  from  that  of  the  corresponding 
response  signals  of  the  heated  wire  to  obtain  the  power  density  spectrum  of  the  “noise-free” 
wire  response  signal. 

The  response  of  each  of  the  principal  instruments  used  to  process  the  a-c  voltage 
component  of  the  hot-wire  anemometer  output  prior  to  spectral  analysis  was  characterized 
by  a  rolloff  with  frequency  within  the  bandwidth  of  interest  for  the  measurements.  The 
method  used  by  Laderman  and  Demetriades  (Ref.  26)  to  enhance  the  output  signals 
mathematically,  to  compensate  for  response  degradation,  has  been  used  in  processing  the 
data  acquired  in  the  present  investigation.  For  this  purpose,  the  rolloff  of  the  response  of  the 
wire  sensor,  without  electronic  compensation,  was  defined  mathematically  as  a  function  of 
frequency.  The  frequency  characteristics  of  the  amplifier  thermal  lag  compensation  circuits, 
as  a  function  of  time-constant  setting,  and  those  of  the  amplifier  high-cutoff  filters  were 
determined  by  calibration  by  the  VKF,  and  the  results  were  expressed  mathematically  for 
each  instrument  as  a  function  of  frequency.  The  inverse  of  each  of  the  mathematical 
expressions  for  (1)  the  response  of  the  wire,  (2)  the  compensation  circuit,  and  (3)  the  filtering 
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circuit  was  applied  as  a  transfer  function  in  the  discrete-point  data  reduction  procedures. 
The  resulting  “restored”  form  of  the  wire  response  spectrum  was  used  in  subsequent  data 
processing. 

Electronic  compensation  of  the  constant-current  anemometer  for  thermal  lag  in  the 
sensing  wire,  using  the  dedicated  stage  of  the  amplifier  (Section  4.3.1),  requires  that  the 
amplifier  be  adjusted  manually  to  the  appropriate  thermal  time  constant  of  the  wire,  a  factor 
that  is  altered  by  a  change  of  the  heating  current  of  the  wire  or  a  change  of  the  unit  Reynolds 
number  of  the  flow  around  the  wire.  (The  appropriate  value  of  the  wire  lime  constant  is 
determined  using  the  square-wave  technique,  Ref.  27.)  However,  appropriate  compensation 
for  sensor  wire  thermal  lag  can  also  be  achieved  mathematically,  through  the  use  of  the 
transfer  function.  This  is  done  after  the  fact  and  in  lieu  of  manual  adjustments  of  the 
amplifier  time  constant,  for  each  heating  current,  before  data  acquisition.  This  option  is 
available  provided  that  the  time  constant  of  the  given  wire  has  been  determined  (by  the 
square-wave  method)  for  one  heating  current  at  one  test  condition.  This  optional  procedure 
was  used  in  the  present  investigation  to  improve  the  efficiency  of  the  measurements  and  to 
minimize  the  exposure  time  of  the  wire  to  the  hazards  of  the  flow  field. 

At  this  point  in  the  signal  processing,  the  power  density  spectrum  of  the  hot-wire  a-c 
voltage  response  would  be  noise  free,  enhanced  to  restore  signal  content  lost  because  of 
instrumentation  deficiencies,  and  compensated  for  wire  thermal  lag.  The  discrete  point 
spectrum  can  at  this  point  be  summed  over  the  frequency  range  to  yield  the  wide-band  mean- 
square  voltage  fluctuation  of  the  wire.  It  was  found  in  the  present  investigation  that  the 
wide-band  fluctuations  of  signals  produced  by  low-level  turbulence  are  well  behaved  when 
processed  in  the  manner  and  sequence  discussed  above.  This  finding  is  in  contrast  with 
results  obtained  using  the  measured  wide-band  signal  levels:  when  the  mean-square  of  the 
measured  wide-band  noise  signal  (unheated  wire)  is  subtracted  from  the  mean-square  of  the 
measured  wide-band  response  signal  of  the  heated  wire,  the  results  are  subject  to 
uncertainties  which  can  become  unacceptable  as  the  level  of  the  response  signal  approaches 
that  of  the  noise. 

From  this  point  in  the  signal  processing,  the  data  reduction  can  proceed  in  two  ways: 
(1)  the  wide-band  (summed  spectrum)  voltage  fluctuation  data  can  be  reduced  to  the  form  of 
wide-band  flow  fluctuations,  and  (2)  the  discrete  point  specira  of  the  voltage  fluctuation 
data  can  be  reduced  frequency  by  frequency  to  the  form  of  discrete  point  spectra  of  the  flow 
fluctuations.  The  modal  analysis  processes  discussed  in  the  next  section  are  applicable  to 
both  the  wide-band  and  the  spectral  data  reduction,  unless  noted  otherwise. 


25 


A  EDC-TR  -81  -5 


5.2.3  Modal  Analysis 

Procedures  used  in  the  interpretation  of  hot-wire  anemometer  measurements  in  a  flow 
field  for  the  purpose  of  inferring  (l)  which  flow  fluctuations  are  present  and 
(2)  the  magnitudes  of  the  fluctuations  are  generally  referred  to  as  “modal  analysis” 
procedures.  In  supersonic  flow  the  hot  wire  responds  to  fluctuations  of  mass  flux  and  total 
temperature  (the  sensing  variables)  whereas  the  fluctuations  which  characterize  an  unsteady 
flow  field  are  those  of  vorticity,  entropy,  and  sound  (the  independent  mode  variables). 
Fluctuations  in  the  sensing  variables  and  the  mode  variables  are  inferred  from  hot-wire  a-c 
voltage  response  (rms)  measurements  using  the  “mode”  concept  and  associated  techniques 
of  “fluctuation  diagrams”  introduced  by  Kovasznay  (Ref.  28)  and  expanded  by  Morkovin 
(Ref.  29).  The  comprehensive  discussion  by  Morkovin  (Ref.  29)  of  the  derivation  of  the 
basic  relations  for  the  response  of  the  hot  wire  has  been  outlined  by  various  authors 
including  Donaldson  and  Wallace  (Ref.  7)  and  has  not  been  repeated  in  the  present  report. 
The  application  of  modal  analysis  procedures  to  hot-wire  measurements  made  in  high-speed 
flow  has  been  discussed  in  detail  by  Laderman  and  Demetriades  (Ref.  26)  with  special 
emphasis  on  flow  fields  characterized  by  the  presence  (significance)  of  more  than  one 
fluctuation  mode  variable  and  the  possibility  of  cross-correlation  between  modes. 
Application  in  the  case  of  a  pure  sound  field,  where  only  the  sound-wave  mode  is  present 
and  can  be  attributed  to  a  single  sound  source,  has  been  discussed  by  Vrebalovich  (Ref.  30). 
The  methods  of  Refs.  26  and  30  have  been  used,  as  applicable,  in  the  present  investigation  to 
accomplish  the  modal  analysis  and  the  subsequent  determination  of  the  fluctuations  of  the 
flow  variables,  namely  velocity,  static  temperature,  pressure,  and  density. 

5. 2.3.1  Wire  Response 

The  mean-square  a-c  voltage  fluctuation  (e)2  of  the  hot-wire  output  may  be  expressed  in 
terms  of  the  three  nondimensional  mode  “amplitudes,”  vorticity  (t),  entropy  (a),  and  sound 
(9),  and  the  three  cross-correlations  between  the  modes  by  the  relation  (Refs.  29  and  7) 

(e)2  =  Ae2(r)2  +•  Ae'2(ff)2  4-  Ae2t9)2  -  2  Ae  Ae  (}) (ct)  N 

r  o  n  t  a  to 

...  .  (13) 

-  2Ae  Ae  4  2Ae  Ac  It 

on  Qlt  TTf  jtt 

where  Aer,  Ae„,  and  Ae,  are  the  vorticity,  entropy,  and  sound  sensitivity  coefficients, 
respectively,  of  the  wire,  where  ( — )  denotes  the  rms  value  of  the  fluctuation,  and  where  each 
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R  is  a  correlation  coefficient  identified  by  its  two  subscripts;  for  example  RT„  =  (t<t)/(t)(5)  is 
the  correlation  coefficient  between  fluctuations  of  vorticity  and  entropy.  In  this  presentation 
factors  affecting  the  rms  a-c  voltage,  such  as  the  d-c  gain,  the  signal  enhancement,  and  the 
wire  time  constant  have  been  accounted  for  in  the  term  e;  moreover,  the  “finite-circuit” 
parameter  (which  accounts  for  the  finite  impedence  of  the  “constant"  current  control)  has 
been  accounted  for  in  each  of  the  sensitivity  coefficients. 

The  sensitivity  coefficients  (Ae(  >)  of  a  given  wire  are  functions  of  its  physical  and  heat- 
transfer  characteristics  (see  Table  1  in  Ref.  7),  which  are  determined  from  the  results  of  the 
oven  and  flow  calibrations  of  the  wire  (Appendix  A).  Changes  of  the  wire  heating  current  or 
local  unit  Reynolds  number  of  the  flow  around  the  wire  produce  changes  in  the  sensitivities, 
so  that  the  values  of  the  coefficients  must  be  determined  at  each  operating  condition.  It 
should  be  noted  that  in  supersonic  flow  the  wire  Nusselt  number  and  temperature  recovery 
factor  are  essentially  independent  of  Mach  number. 

In  principal,  six  equations  in  six  unknowns  (the  three  mode  “amplitudes”  and  the  three 
correlation  coefficients)  could  be  obtained  for  a  simultaneous  solution  if  measurements  were 
made  at  six  distinct  wire  heating  currents.  In  reality,  the  sensitivity  coefficients  and  the  rms 
voltage “e  have  uncertainties  which  rule  out  this  method  of  solution.  It  must  be  noted  that  the 
presence  of  sound-wave  fluctuations  involves  amplitude  and  direction  of  propagation  (that 
is,  two  possible  variables,  rather  than  one).  For  simplicity  in  the  present  discussion  it  has 
been  assumed  that  there  is  a  preferred  direction  of  propagation  of  the  sound  waves  and  that 
the  direction  is  independently  known,  unless  otherwise  indicated. 


5.2.3. 2  Fluctuation  Diagram  Technique 


The  fluctuation  diagram  technique  makes  use  of  the  fact  that  the  relative  sensitivity  oT 
the  hot  wire  to  the  various  flow  parameters  varies  with  the  temperature  of  the  heated  wire 
relative  to  that  of  the  unheated  wire  (that  is,  with  wire  “overheat").  By  making 
measurements  at  different  wire  overheat  temperatures  one  obtains  a  diagram  of  the  wire 
response  which  has  a  characteristic  shape  for  each  mode  variable.  More  generally,  when 
modes  coexist  and  arc  cross-correlated,  the  diagram  technique  can  be  used  —  but  only  under 
certain  restrictive  conditions. 


For  supersonic  flow  the  relation  Eq.  (13)  can  be  rewritten  as  follows: 

V2'i2-r2X2  ,  *2(g|  tg2X^-2?iRrffX+2(g|..  RgX)  -  *rRffrX>  (I4> 
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where 


and 

with 


and 


V  =  c/Ae 

a 

AeW  /  Aerr  =  «]  +  VX 

-ty  -  L)  ( iV] 2  -  1) 

8.  ‘ - 

1  i-  (y  -  l)M2 

y(l  +  nxM)  -  v\!2  -  D  tn^/Mi 
1  +  (y  -  1)  VI2 


(15) 

(16) 

(17) 

(18) 

(19) 


where  nx  is  the  direction  cosine  of  the  normal  to  the  sound  plane  wave  front  relative  to  the 
local  mean  flow  direction.  (The  hot  wire  aligned  normal  to  the  mean  flow  direction  is 
sensitive  to  velocity  variations  in  that  direction,  only.) 


Equation  (14)  can  be  expanded  in  terms  of  X  to  yield 


Y2-(52  +  gJff2  +  2giairllffff)T2lr'illrff 

1  6,^r'[^rH  g,?2^2)X  +(rz  +  g2 n2  +  2g2^ll7rr)X2 


(20) 


which  is  of  the  form 

Y2  =0  i  D  X  +  15  X2  (2D 

1  2  3 

where  the  coefficients  D,  represent  various  combinations  of  the  three  mode  amplitudes  and 
their  cross-correlations.  If  one  mode  variable  predominates  in  the  expressions  for  the 
coefficients  so  that  the  remaining  terms  are  negligible,  then  the  form  of  Eq.  (20)  degenerates 
to  a  pair  of  straight  lines,  as  follows: 

Y  =  ±(13  h-  15.  X)  (22) 

r  4  5 

Operation  of  the  hot  wire  at  a  fixed  location  in  the  flow  field  for  several  wire  heating 
currents  (up  to  1 1  in  the  present  work)  will  yield  corresponding  values  of  the  sensitivity 
coefficients  and  of  the  rms  voltage,  e.  For  each  heating  current  one  point  in  the  X-Y  plane 
can  be  determined  using  Eqs.  (15)  and  (16).  The  set  of  points  (X,  Y)  for  the  several  heating 
currents,  plotted  in  rectangular  coordinates,  will  indicate  whether  the  fairing  is  of  the  form 
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of  Eq.  (21)  or  Eq.  (22).  Values  of  the  coefficients  of  the  appropriate  fairing,  obtained  using 
the  method  of  Least-squares,  are  matched  with  the  corresponding  expressions  of  the 
coefficients. 

5. 2.3. 3  Single  Mode 


The  following  discussion  will  illustrate  the  case  of  one-mode  predominance: 

(1)  If  vorticity  (turbulent  velocity)  fluctuations  predominate  so  that  o  =*  ?r  =  0,  Eq .  (20) 
becomes 


Y^x 


(23) 


(2)  If  entropy  (temperature  spottiness)  fluctuations  predominate  so  that  r  *  t  =s  0, 
Eq.  (20)  becomes 

Y  -  to  (24) 

(3)  If  sound  (sound  wave)  fluctuations  predominate  so  that  7  =  a  *  0,  Eq.  (20)  becomes 


Y-±<e,Hg2X);  (25) 

If  the  “best”  fairing  among  a  set  of  points  (X,  Y)  is  a  straight  line  with  nonzero  intercept 
and  nonzero  slope,  the  fluctuation  diagram  technique  would  indicate  the  wire  was  located  in 
a  pure  sound  field.  The  principal  variable  in  this  case  is  pressure,  and  the  relations  among 
the  thermodynamic  variables  are  isentropic.  If  the  least-squares  fairing  is  given  by 


Y  =  13  +  D  ,X 

4  5 

and  the  coefficients  are  matched  with  those  of  Eq.  (25),  the  amplitude  of  the  sound  mode  is 
found  by 

"=UYV  (26) 


From  Morkovin  (Ref.  29)  with  a  =>  0,  fluctuation  variables  can  be  determined  using  the 
relations 


-  =  p  _  P  =  i  T 
p  y  p  >'  -  1  T 


(27) 


A  matching  of  the  slope  coefficients  yields 


(28) 
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If  Eq.  (19)  is  solved  for  the  direction  cosine  nx,  one  obtains  in  the  test  section  flow 


nv  =  bA- 


(v  -  1)  Ml  1-  I 


(29) 


If  g2  is  obtained  from  Eq.  (28)  and  substituted  into  Eq.  (29),  a  preferred  orientation  of  the 
sound  waves  intercepted  by  the  hot  wire  can  be  determined.  Laufer  (Ref.  31)  has  adapted  the 
theory  of  Phillips  (Ref.  32)  for  the  radiation  of  sound  by  a  free  shear  layer  to  the  case  of 
boundary- layer  sound  radiation.  In  essence,  a  frozen  eddy  pattern  acting  as  a  wavy  wall 
moves  downstream  within  the  boundary  layer  at  a  supersonic  speed  relative  to  the  free- 
stream  flow  and  radiates  energy  in  the  form  of  Mach  waves.  If  the  mean  propagation 
velocity  (streamwise)  of  the  sound  source  is  designated  by  us,  the  mean  velocity  of  the  free 
stream  relative  to  the  virtual  wavy  wall  is  ur  =  u«  -  us,  and  the  corresponding  Mach 
number  is  Mf: 


S.2.3.4  Coexisting  Modes 


If  the  “best”  fairing  among  a  set  of  points  (X,  Y)  is  a  hyperbola,  of  the  form  of  Eq.  (21), 
a  matching  of  the  coefficients  with  those  of  Eq.  (20)  will  give  a  set  of  three  equations  in  the 
seven  unknowns  (including  directior)  of  sound  wave  propagation).  Resolution  of  the 
fluctuations  of  coexisting  modes  from  this  underdetermined  set  of  equations  requires  more 
information  from  independent  sources,  for  example;  (1)  an  acceptable  fluid  dynamics  model 
of  the  flow  field  around  the  wire  which  would  supply  the  required  number  of  additional 
independent  relations  in  some  (or  all)  of  the  unknowns  and/or  (2)  measurements  of  any  of 
the  fluctuations  by  independent  methods  which  would  serve  to  relax  the  a  priori  constraints 
inherent  in  use  of  the  fluid  dynamics  model.  In  the  present  discussion  three  models  suggested 
by  Laderman  and  Demetriades  (Ref.  26)  have  been  considered: 

Model  1:  No  sound;  that  is,  nr  =  0 


with  this  constraint,  Eq.  (20)  reduces  to 

V"  -  \  2rr;H.f;\  -  t 

and  matching  coefficients  with  the  fitted  hyperbola  will  yield  values  for  vorticity  and 
entropy  fluctuations  and  for  their  correlation  coefficient.  This  model  is  often  used  in 
reducing  hot-wire  anemometer  measurements  obtained  in  supersonic  turbulent  boundary- 
layer  How  where  sound  fluctuations  are  considered  negligible  compared  to  those  of  vorticity 
and  entropy. 
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Model  2:  Uncorrelaied  modes;  that  is, 


With  these  constraints,  Eq.  (20)  reduces  to 

Y 2  =  (  »2  +  g2 ;2)  +  2g]g^2X  +  ( r'2  +  ^2)\2  (32) 

A  matching  of  coefficients  with  the  fitted  hyperbola  can  be  made  to  yield  a  completely 
determined  set  of  equations,  if  a  preferred  orientation  of  the  sound  waves  [which  enter 
through  the  direction  cosine  nx  in  the  term  g2,  as  shown  in  Eq.  (19)]  can  be  determined 
independently  (or  assumed).  Kovasznay  (Ref.  28)  and  Morkovin  (Ref.  29)  suggest  that  this 
model  may  be  applicable  in  the  free  stream  of  supersonic  wind  tunnels. 

Model  3:  Uncorrelated  sound  and  anticorrelated  vorticity  and  entropy;  that  is, 

n  -n  -0  and  it  =-1.0 

<  JTT  ? T  T  Ti 7 

With  these  constraints,  Eq.  (20)  reduces  to 


Y2  *  (v2  +  g2^2)  -  X-W+  g1g2?j2)X  +  (f2  +  g2!72)X- 


(33) 


A  matching  of  coefficients  with  the  fitted  hyperbola  can  be  made  to  yield  a  completely 
determined  set  of  equations  if  the  direction  cosine  nx  is  supplied. 


5.3  DYNAMIC  PRESSURE  ANALYSIS 


in  general,  the  PSD  values  presented  herein  are  in  engineering  units,  namely  (psi)2/Hz  or 
(mv)2Hz,  or  are  normalized  as  follows: 

V/^lu.,/8*  ■  (34) 

Note  that  C(  )  is  simply  a  pressure  coefficient  formed  as  a  ratio  of  the  fluctuating  pressure 
parameter  relative  to  the  mean  pressure  level  divided  by  the  free-stream  dynamic  pressure. 

A  Sirouha!  number  is  formed  for  the  frequency  coordinate  of  the  spectral  distributions, 
as  follows: 


Strouhal  Number  =  f  •  6*/u» 

If  the  convective  velocity  uc  of  the  sound  source  associated  with  the  tunnel  wall  turbulent 
boundary  layer  had  been  well  established,  it  would  have  been  a  more  acceptable  parameter 
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to  use  with  the  aforementioned  parameters  than  the  free-stream  velocity  component  (see 
Section  2.0). 

The  remaining  major  normalizing  parameters  concern  the  ratio  of  the  accumulated 
energy  in  the  spectrum  up  to  a  given  frequency  to  the  total  energy  in  the  spectrum.  This 
normalized  energy  distribution  (NED)  is  defined  as  follows: 


k 

N£D(k  •  .\D  -  £ 

n=L 


f  a/2 

Hu • AO /  £ 

/  n=l 


Rn  •  AO 


Note  that  BW  =  (N/2)  •  Al\ 


(35) 


The  frequency  at  which  this  ratio  reaches  (for  example,  0.95)  simply  identifies  the  signal 
bandwidth  of  the  spectrum  which  contains  95  percent  of  the  energy  of  the  fluctuation 
parameter. 

Two  fundamental  relationships  are  used  in  this  analysis  in  an  attempt  to  correlate  the 
results  of  the  hoi-wirc  aneinometry  to  the  dynamic  pressure  measurements.  As  noted  in  Ref. 
14,  for  a  plane  wave,  the  free-stream  pressure  and  velocity  fluctuations  are  related  as 
follows: 

P/P  =  yM^tii/u)  (36) 

or 

C-  =  2tu/Tf) 

The  pressure  measurements  obtained  from  a  pitot  probe  instrumented  with  dynamic 
pressure  transducers  is  related  to  the  sLatic  pressure  fluctuations  deduced  from  the  hot-wire 
probe  (Ref.  .33).  A  good  approximation  to  the  relationship  of  the  pitot  probe  pressure  and 
Lhe  free-stream  properties  is  as  Follows: 


where  the  constant  “A”  is  primarily  a  function  of  y,  the  ratio  of  the  specific  heats  of  the 
gas,  and  a  very  weak  function  of  the  free-stream  Mach  number.  Actually,  the  value  “A” 
varies  from  0.953  at  Mach  number  1 .5  to  1.085  at  Mach  number  8  (y  =  1.4)  and  hereafter  is 
arbitrarily  set  equal  to  one. 


Equation  (37)  can  be  expanded  in  terms  of  the  instantaneous  fluctuations  relative  to  the 
mean  value,  neglecting  higher  order  terms,  as  follows: 


APp/P,  -  (_V /p)  +  2(j_\li /ul 


(38) 
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Forming  the  mean  square  of  this  relationship  of  the  instantaneous  values  yields  the 
following  result: 


(APp/Pp)2  "  (iV/p)2  +  l('Wu>2  +  -KAp/pKAu/u) 

Typically, 

S!wr'pJ*il 

where  (Pp/Pp)2  is  the  ratio  of  the  mean  square  value  of  the  fluctuation  relative  to  the  mean 
value,  P p.  But  the  last  term  on  the  right  in  Eq.  (39)  is  the  cross  correlation  between  the  free 
stream  density  and  velocity  fluctuations  and  may  be  positive  or  negative.  In  general  the 
averaged  (mean)  value  of  this  term  is  defined  as  follows: 


As  noted  in  Ref.  14  the  Correlation  factor  was  identified  as  being  negative  since  the 
fluctuating  properties  are  assumed  to  behave  isentropically  and  the  disturbance  is  primarily 
the  result  of  a  pure  sound  Field  emanating  from  the  tunnel  wall  boundary  layer.  In  addition, 
the  fluctuations  can  be  viewed  as  a  plane  wave  with  the  following  characteristics: 


Thai  is, 


and,  as  noted  in  Ref.  14, 


(41) 


nx  =  l[(u8/u)  -  1MJ-1  (42) 

Therefore, 

■Jiyiy2-  2  u  -  2n(/wj  («> 

Based  on  the  present  hot-wire  analyses  and  previous  experimental  work  (Ref.  14),  the 
ratio  of  the  source  velocity  to  the  free-stream  velocity  was  set  equal  to  a  nominal  value  of 
0.40,  Therefore,  Eq.  (43)  was  simplified  as  follows: 
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As  noted  in  Ref.  34  and  included  in  the  previous  expression,  the  measured  fluctuating  pitot 
probe  pressure  is  twice  the  actual  value  in  order  to  account  for  the  reflection  of  the 
impinging  fluctuating  pressure  wave  on  the  probe  face. 

The  terms  APp/P2,  A q/q2,  and  Au/u2  in  Eq.  (39)  can  be  evaluated  as  the  mean  square 
value  of  the  deviations  from  the  mean  value  of  the  parameter,  but  the  term  (Ae/p)(Au/u)  is 
a  cross-correlation  term  which  may  be  either  positive  or  negative.  As  noted  in  Ref.  14  and  as 
exhibited  in  the  modal  analysis  of  the  test  data,  this  term  is  negative  if  the  disturbance  is 
assumed  to  be  isentropic  and  the  result  of  a  pure  sound  field.  In  addition,  the  fluctuations 
can  be  viewed  as  a  plane  wave  with  the  following  characteristics: 


An  _  %  Ap 

yp  (45) 


Basically,  this  analysis  assumed  that  the  rms  value  defined  in  Eq.  (41)  can  be  redefined 
simply  as 


(46) 


i 

In  this  form  the  magnitude  and  the  sign  of  the  direction  cosine  (nx)  between  (u/u) 
and  p/p)  are  preserved  and  unaffected  by  the  rms  procedure.  If,  on  the  other  hand,  the 
absolute  value  of  nx  must  prevail,  the  agreement  between  the  hot-wire  and  the  pitot  acoustic 
results  becomes  extremely  poor  at  the  lower  free-stream  Mach  numbers  since  Eq.  (44)  must 
be  defined  as  follows: 


Correspondingly,  Eq.  (43)  becomes 
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6.0  RESULTS 


6.1  TUNNEL  A  STILLING  CHAMBER 

The  stilling  chamber  measurements  in  Tunnel  A  were  made  downstream  of  the  screens 
near  the  maximum  diameter  of  the  stilling  chamber  (see  Fig.  Id).  The  pressure  fluctuation 
levels,  as  expected,  are  essentially  independent  of  the  nozzle  throat  contour  but  appear  to  be 
directly  related  to  the  mean  value  of  the  total  pressure  in  the  stilling  chamber,  as  shown  in 
Fig.  15.  The  addition  of  high-pressure  makeup  air  upstream  of  the  screens  (using  valve  337 
or  237)  momentarily  perturbs  this  level  to  a  higher  value  by  more  than  an  order  of  magnitude 
(an  increase  of  nominally  24  db  for  the  case  shown  in  Fig.  15).  Except  for  the  effects 
produced  by  the  make-up  air,  the  stilling  chamber  pressure  fluctuation  levels  appear  to  vary 
with  the  measured  stilling  chamber  pressure  (PT)  as  follows: 

PT/PT  (percent)  =  0.235(PT)  90  (49) 

for  4  psia  >  PT  >  130  psia  in  Tunnel  A. 

A  comparison  of  some  previously  unpublished  hot-wire  results  obtained  in  November 
1965  with  the  present  dynamic  pressure  measurements  obtained  in  Tunnel  A  stilling  chamber 
is  shown  in  Fig.  16.  The  temperature  fluctuation  levels  based  on  the  hot  wire  anemometer 
measurements  were  negligible,  less  than  0.1  percent  of  the  stilling  chamber  mean  value.  In 
Ref.  4,  the  relation  between  the  velocity  fluctuations  obtained  from  the  hot-wire 
measurements  and  those  deduced  from  a  stilling  chamber  surface-mounted  dynamic 
pressure  gage  is  as  follows: 

'W  ■  <50> 

Using  this  relationship  and  converting  the  measured  stilling  chamber  pressure  fluctuation  to 
the  corresponding  velocity  values  shows  (Fig.  16b)  that  the  hot-wire  results  obtained  in 
November  1965  (open  symbols)  are  in  fair  agreement  with  the  present  pressure  fluctuation 
measurement  (solid  symbols).  The  flagged  values  correspond  to  fluctuation  levels  occurring 
at  the  maximum  tunnel  stagnation  pressure.  In  general,  the  velocity  fluctuations  were  found 
to  be  less  than  1  percent,  and  the  measured  level  of  the  pressure  fluctuations  was  nominally 
less  than  0.1  percent. 
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Samples  of  the  spectra  of  the  pressure  fluctuation  signal  obtained  in  the  stilling  chamber 
are  given  in  Fig.  17  for  various  stilling  chamber  pressure  levels  and  plant  staging 
configurations.  Also,  the  frequencies  where  more  significant  energy  levels  occur  are 
identified.  Most  of  the  energy  occurs  below  5  kHz.  A  comparison  of  Figures  17e  and  f 
clearly  shows  that  the  addition  of  makeup  air  to  maintain  the  tunnel  stilling  chamber 
pressure  level  (flow  rate  or  test  section  free-stream  Reynolds  number)  produces  a  significant 
disturbance  which  in  previous  acoustical  test  programs  has  compromised  the  dynamic 
pressure  signals  recorded  in  the  tunnel  test  section. 

6.2  TUNNEL  A  AND  B  TEST  SECTION  TURBULENCE 

The  evaluation  or  the  test  section  turbulence  levels  in  Tunnels  A  and  B  is  based  on  free- 
stream  hot-wire  anemomelry  measurements,  free-stream  dynamic  pitot  pressure  probe  data 
(pitot-acoustic  probe),  and  surface  static  pressure  fluctuation  measurements.  The  surface 
pressure  fluctuation  measurements  were  obtained  from  a  Hat  plate  and  the  Tunnel  A  side 
wall  surfaces  in  1978  and  from  a  6-deg  sharp  cone  in  1979. 

The  1978  flat  plate  surface  measurements  were  subsequently  found  to  be  unsatisfactory 
(signal-to-noise  ratio  was  too  low}  for  defining  the  test  section  turbulence.  The  nylon  insert 
containing  the  dynamic  transducers  (see  Fig.  9)  was  too  flexible  and,  when  pressure  was 
applied  to  the  nylon  surface,  tended  to  preload  or  alter  the  signal  level  from  the  piezoelectric 
dynamic  pressure  sensor.  The  signal  from  the  smaller  transducer  was  extremely  noisy  and 
indicated  that  the  gage  was  not  flush  with  the  nylon  insert  (i.  e.,  the  mid-frequency  range  of 
the  signal  was  attenuated).  Also,  the  Corcos  gage  size  correction  factor  (sec  Ref.  35)  for  the 
two  transducers,  nominally  0.2  and  0.4  in.  in  diameter,  was  extremely  large.  It  seemed  to  be 
unrealistic  to  apply  these  gage  size  corrections  to  such  poor  signal  responses. 

In  all  subsequent  installations,  the  dynamic  pressure  transducers  were  mounted  in  a  thin 
cylindrical  RTV  sleeve  to  insulate  the  gages  from  vibrational  (or  acceleration)  disturbances 
and  to  isolate  potential  acoustical  noise  sources  to  prevent  their  being  transmitted  through 
the  model  to  the  gage.  Although  the  present  flat  plate  results  are  highly  questionable  and  are 
not  included  in  this  document,  this  flat  plate  acoustics  model  is  available  for  future 
investigations. 

6.2,1  Test  Section  Pressure  Fluctuations 

A  comparison  of  the  free-stream  static  pressure  fluctuations  inferred  from  the  hot-wire 
anemomelry  and  pilot-acoustic  probe  data  as  a  function  of  frcc-stream  Reynolds  number 
based  on  the  effective  test  section  diameter  (DE)  is  given  in  Figs.  18  and  19.  In  general, 
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considering  the  data  scatter,  the  resuhs  from  the  two  measuring  techniques  are  in  agreement 
in  Tunnel  A  (Fig.  18)  and  in  Tunnel  B  at  the  maximum  Reynolds  number  (Fig.  19).  The 
difference  in  the  results  between  the  two  measurement  techniques  at  Mach  number  6  at  the 
low  Reynolds  numbers  is,  at  this  time,  unresolved. 

In  general,  the  results  obtained  at  Mach  numbers  2.0,  4.0  and  5.0  in  Tunnel  A  and  the 
pitot  acoustic  data  in  Tunnel  B  indicate  that  the  static  pressure  fluctuations  vary  as 
(RE|3E)^i/4,  which  is  in  reasonable  agreement  with  the  data  trends  shown  in  Ref.  5.  The  hot 
wire  resuhs  in  Tunnel  B  do  not  follow  this  Reynolds  number  trend  (that  is,  Cp  = 
f(REDE_  1/4).  Although  the  actual  reason  for  this  anomaly  is  presently  unknown,  it  is 
possible  that  the  difference  in  the  hot-wire  and  pitot  acoustic  probe  location  or  orientation 
with  respect  to  the  tunnel  flow  or  tunnel  walls  affected  these  results. 

6.2. 1.1  Hot-Wire  Anemometer  Spectral  Results 

Selected  spectra  of  the  fluctuations  determined  from  the  hot-wire  measurements  are 
presented  in  Figs.  20  through  24.  These  spectra  were  obtained  with  a  free-stream  unit 
Reynolds  number  of  RE*  =  3  million  per  foot  for  M*  =  4,  6,  and  8. 

At  M„  =  4  and  6  the  fluctuations  in  the  free-stream  flow  are  dominated  by  the  sound 
mode  (Section  5. 2. 3.4),  wherein  the  principal  variable  is  pressure.  Velocity  and  static 
temperature  fluctuations  in  this  case  are  induced  by  the  pressure  fluctuations  and  are 
inferred  from  them  using  isentropic  relations,  which  are  simply  multiplication  factors.  (The 
sound  waves  are  assumed  to  have  a  single  or  “preferred”  direction.)  As  a  consequence,  the 
pressure  fluctuation  spectrum  serves  to  illustrate  the  spectra  of  the  other  fluctuating 
parameters. -The  speci rum  of  the  pressure  fluctuations  for  RE*  =  3  million  per  foot  at  M» 
=  6  is  shown  plotted  in  a  semilog  format  in  Fig.  20,  and  the  spectrum  for  the  same  unit 
Reynolds  number  at  M*  =  4  is  shown  in  a  log-log  format  in  Fig.  21. 

The  power  density  spectrum  of  the  pressure  fluctuations  for  M*.  =  4  is  presented  in 
Fig.  22  with  the  amplitude  normalized  using  the  value  of  P/P  =  0.000669  (for  f  =  150  Hz), 
which  is  the  value  of  P/P  at  150  Hz.The  frequency  has  been  nondimensionalized  using  the 
free-stream  velocity  and  the  longitudinal  Integral  scale  given  by 

L  P(f  =  150  El  ) 

A  =  —  — - — 

P  *  P(\s  i  rl o  biind) 

The  velocity  fluctuation  power  density  spectrum  normalized  using  u/u  =  0.0000445  (for 
f  =  150  Hz)  is  also  shown  in  Fig.  22,  where  the  integral  scale 
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has  been  used  in  nondimensionalizing  the  frequency. 

The  corresponding  spectra  obtained  from  the  M»  =  6  measurements  are  presented  in 
Fig.  23.  In  this  case  the  values  P/P  =  0.00306  and  u/u  =  0.00159  (for  f  =  600  Hz,  in  both 
cases)  were  used  to  normalize  (he  respective  amplitudes,  and  the  integral  scales 

u  p(  f =f.U0  H  ) 

^  _  OC  7. 

P(.  ide  band) 

and  u  u(  f-  600  II  ) 

A  -  - — 

U  ^  u(wide  band) 

were  used  in  nondimensionalizing  the  frequency  in  the  respective  abscissas. 


At  Mo,  =  8  the  free-stream  fluctuations  are  not  dominated  by  the  sound  mode 
(Section  5. 2. 3.4).  Presented  in  Fig.  24  are  the  power  density  spectra  of  the  pressure,  velocity, 
and  static  temperature  fluctuations  obtained  after  it  was  specified  that  1)  the  sound  is 
uncorrelated  with  either  vorticity  or  entropy,  and  that  the  latter  two  modes  are  perfectly 
anticorrelated  (that  is,  R„  =  -  1 .0,  =  Rxr  =  0),  and  2)  the  sound-wave  source  velocity 

is  the  same  at  both  =  6  and  8.  The  fluctuations  have  been  normalized  using  the 
following  values:  P/P  =  0.00452,  u/u  =  0.000164,  and  T/T  =  0.00388  (all  for  f  =  600 
Hz).  The  frequency  has  been  made  nondimensional  using  the  free-stream  velocity  and  the 
respective  integral  scales 
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6.2. 1.2  Pitot- Acoustic  Probe  Spectral  Results 


Typical  power  spectral  density  (PSD)  of  the  pitot  pressure  fluctuations  in  Tunnels  A  and 
B  are  shown  in  Fig.  25.  The  local  free-stream  velocity  and  dynamic  pressure  and  the  tunnel 
wall  boundary-layer  displacement  thickness  (defined  in  Appendix  C)  were  used  to  normalize 
these  spectra.  The  free-stream  Reynolds  number  for  the  spectra  shown  for  Mach  number  1.5 
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(Tunnel  A)  to  6.0  (Tunnel  B)  vary  nominally  from  3  x  106  to  4  x  I06  per  foot.  The 
Mach  number  8  data  of  Tunnel  B  were  obtained  at  a  Reynolds  number  of  1  x  10*  per  foot 
(the  analysis  for  the  higher  Reynolds  number  results  was  questionable).  Except  for  the 
Mach  number  8  results,  the  other  spectra  fall  nominally  within  an  amplitude  bandwidth  of 
±6.5  db. 

The  results  in  Fig.  26a  include  a  repeat  of  the  signal  analysis  obtained  for  the  1979  probe. 
The  difference  between  these  two  spectra  is  nominally  3  db  or  slightly  more  than  2  standard 
deviations  ot  the  random  error  cited  in  Section  5.1.2.  The  effect  of  small  changes  in  free- 
stream  Reynolds  number  on  the  spectra  are  shown  in  Figs.  26b  and  c  for  Mach  numbers  2 
and  3,  respectively.  In  general,  the  Reynolds  number  effects  arc  not  clearly  evident  because 
of  the  statistical  uncertainty  in  the  spectral  analysis.  The  Mach  number  5  spectral 
comparisons  of  the  two  probes  were  deleted  from  this  set  because  the  1978  results  were 
questionable.  Figure  27  contains  a  spectrum  obtained  for  the  1979  pitot  acoustic  probe  in 
Tunnel  B  at  Mach  number  6. 

A  comparison  of  (he  pitot  pressure  spectra  obtained  with  two  different  pitot  acoustic 
probes  (the  1978  and  1979  probes  shown  in  Fig.  13)  is  also  shown  in  Fig.  26.  In  general,  the 
output  from  the  1978  probe  was  highly  susceptible  to  various  high-frequency  vibrations  of 
the  probe  support  system  as  indicated  in  Figs.  28  and  29. 

Although  the  true  rms  (TRMS)  of  the  signal  shown  in  Fig.  28  from  the  FY78  probe  in 
Tunnel  A  was  only  0.1  percent  of  the  atmospheric  pressure,  there  is  a  pronounced  peak  or 
resonance  in  the  spectra  at  320  to  400  Hz  obtained  with  the  probe  in  the  model  injection  tank 
and  with  no  flow  through  the  tunnel  circuit.  This  result  was  unaffected  by  whether  or  not  the 
tunnel  was  operating.  Venting  the  tank  to  the  tunnel  static  conditions  or  reducing  the  model 
injection  tank  pressure  doubled  the  signal  level  and  revealed  additional  vibrational  modes 
which  were  picked  up  by  the  pitot  acoustic  transducer.  Theoretically,  the  harmonics  for 
these  vibrational  bending  modes  vary  as  follows:  6.26  f0,  17.55  f0,  34.39  fG,  etc.  (see  Ref.  36, 
pp.  71-73).  In  general,  the  fundamental  mode  (f0)  is  proportional  to  the  radius  of  gyration 
of  the  beam  or  probe  section  and  inversely  proportional  to  the  square  of  the  beam  (probe) 
length. 

The  results  in  Fig.  28  were  used  as  a  tare  run  to  identify  those  signal  frequency  modes 
which  are  no(  considered  to  be  a  part  or  ihe  signal  when  the  probe  is  exposed  to  the  tunnel 
flow.  As  shown  in  Fig.  26,  the  spectra  were  faired  through  at  these  frequencies  which  were 
attributed  to  transducer  vibration.  The  signals  were  then  reintegrated  to  define  the  true  Fms 
(TRMS)  of  Ihe  pressure  fluctuations  sensed  by  the  probe  when  exposed  to  the  tunnel  flow'. 
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The  pitot  acoustic  probe  installation  in  Tunnel  B  differed  from  that  in  Tunnel  A,  and 
this  resulted  in  different  probe  vibrational  bending  inodes,  as  shown  in  Fig.  29.  The  probe 
consists  of  two  cylindrical  sections  pinned  together.  Apparently  one  clearly  identifiable 
vibration  was  produced  by  the  entire  probe  support  “L2”  and  another  was  produced  by  the 
section  closest  to  the  probe  support,  “LI.”  The  longer  beam  length  was  responsible  for  the 
lower  frequencies. 

As  indicated  in  Fig.  29,  with  no  flow  in  the  circuit  and  the  model  injection  tank  at 
atmosphere,  the  two  fundamental  vibrational  modes  appear  to  be  present  even  though  the 
overall  true  rms  energy  level  signal  is  negligible  (TRMS  <  0.05  percent).  With  flow  in  the 
tunnel,  the  energy  level  (TRMS)  increases  to  0.6  percent.  Venting  the  model  injection  tank  to 
the  tunnel  (reducing  the  tank  pressure)  significantly  increases  the  signal  level  so  that  the 
higher  vibrational  frequencies  appear  in  the  signal  spectrum.  As  in  the  case  of  the  Tunnel  A 
results,  these  vibrational  frequencies  were  then  used  to  fair  through  or  filter  the  spectra 
obtained  when  the  probe  was  exposed  to  the  tunnel  flow. 

6.2.2  Tunnel  A  Side  Wall  Static  Pressure  Fluctuations 

In  Tunnel  A,  two  side  wall  static  taps  were  instrumented  with  dynamic  pressure  sensors. 
The  results  shown  in  Fig.  30  identify  the  total  amount  of  energy  existing  (NED)  in  the  profile 
up  to  a  given  frequency.  For  example,  of  the  energy  detected  by  the  nonoperating  side  wall 
(the  side  opposite  the  tunnel  control  room)  gage,  80  percent  is  contained  in  the  frequency 
envelope  up  to  22.5  kHz.  For  the  operating  side,  80  percent  of  the  energy  detected  is  in  the 
frequency  envelope  up  to  12.5  kHz. 

Although  the  two  gages  are  supposedly  identical,  the  operating  side  gage  must  be 
protruding  slightly  from  the  tunnel  wall,  causing  a  reduction  in  the  mid  range  frequency 
response  of  the  gage  as  shown  in  Fig.  30b  (see  Ref.  17).  Unfortunately,  the  gage  installation 
was  removed  before  this  problem  with  the  operating  side  wall  gage  could  be  visually 
confirmed.  Obviously,  the  capability  of  performing  some  online  spectral  analysis  would  aid 
in  the  validation  of  dynamic  measurements. 

The  wall  static  spectra  obtained  on  the  nonoperating  tunnel  side  wall,  shown  in  Fig.  31 , 
provide  a  reasonable  description  of  the  eflect  of  Mach  number  on  the  tunnel  wall  static 
pressure  fluctuations  in  Tunnel  A.  The  levels  of  these  speclra  are  higher  than  those 
associated  with  the  free-stream  lest  section  pressure  fluctuation  levels  deduced  from  the  hot¬ 
wire  or  pitot  acoustic  probes.  These  results  in  Fig.  31  are  simply  a  fairing  of  the  spectra 
obtained  at  various  Tunnel  A  Mach  number  nozzle  settings  at  a  nominal  Reynolds  number 
of  4  x  106  per  foot. 
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6.3  CONE  SURFACE  STATIC  PRESSURE  FLUCTUATIONS 

The  6-deg  sharp  cone  was  instrumented  with  dynamic  pressure  gages  at  nominally  three 
stations  along  the  model.  The  rms  level  is  shown  in  Fig.  32  in  coefficient  form  as  a  function 
of  free-stream  Reynolds  number  based  on  the  axial  distance  from  the  nose  to  the 
instrumented  station  (Sta).  These  results  can  be  used  to  identify  the  nominal  level  of  the 
surface  pressure  fluctuation  sensed  on  a  model  in  Tunnels  A  and  B  under  a  laminar, 
transitional,  and  in  some  cases,  a  fully  turbulent  boundary  layer. 

In  general,  as  shown  in  Fig.  32,  the  larger  size  gages  (square  symbols)  produce  an 
attenuated  response  relative  to  the  smaller  gages  (circular  symbols).  This  reduction  in  the 
response  of  the  larger  gage  is  expected  and  simply  reflects  that  a  correction  like  the  one 
formulated  by  Corcos  (Ref.  35)  and  discussed  in  Appendix  B  must  be  applied  before  the  true 
pressure  fluctuation  level  is  experimentally  established.  The  applicability  of  the  Corcos  type 
of  correction  to  results  obtained  in  the  environment  of  a  large  supersonic/hypersonic  test 
unit  would  have  to  be  thoroughly  evaluated  before  it  could  be  applied  as  a  standard 
procedure.  Presently,  the  minimum  size  transducers  used  in  Tunnels  A  and  B  are  0. 10  in.  in 
diameter,  which  means  that  the  maximum  ratio  of  (u  •  r/ut)  varies  from  1.5  to  1.8  (for  f  = 
100  kHz)  depending  on  the  free-stream  Mach  number  (and  the  corresponding  convective 
velocity).  Large  values  of  oj  *r/uc  result  in  significant  gage  attenuation  effects. 

Except  for  the  results  at  Mach  number  2,  these  results  indicate  that  the  nominal 
transitional  Reynolds  number  (where  Cp  reaches  a  maximum  value)  is  nominally  4  x  106  and 
is  essentially  independent  of  the  free-stream  Mach  number.  In  general,  the  transition 
location  based  on  past  observations  of  a  sharp  cone  (half-angles  less  than  9  deg)  can  be 
empirically  correlated  as  follows  for  Tunnels  A,  B  and  C  with  the  model  wall  temperature 
nominally  at  540  to  600°R: 

X|(end)  =  < . 5(RI',/in.  ■  10“ 4  (51a) 

X((. beginning)  =  4.6<HE/in.  •  lO"6)^  (51b) 

This  correlation  can  be  rearranged  and  solved  to  show  that  the  transition  locations  indicated 
in  Fig.  32  arc  bracketed  by  the  correlation. 

In  general,  the  gages  of  the  same  size  located  diametrically  opposite  each  other  at  model 
station  1 1 .0  produced  the  same  result.  Small  differences  between  the  two  measurements  can 
be  attributed  to  slight  deviations  in  the  pitch  angle  from  the  desired  zero  angle  of  attack 
attitude,  which  has  an  effect  on  the  state  of  the  boundary  layer.  Also,  uncertainties 
associated  with  the  gage  calibration  and  flushness  of  the  gage  installation  can  affect  the 
transducer  response  and  contribute  to  the  differences  shown. 
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Figure  33  illustrates  the  effects  of  free-stream  Reynolds  number  (boundary-layer 
transition)  and  gage  location  on  the  power  spectral  density  distribution  of  the  pressure 
fluctuation  measurements  [(p$i)2/Hz  versus  frequency]  made  on  the  6-deg  sharp  nose  cone 
at  Mach  number  6.  An  interesting  phenomenon  appeared  in  the  spectra  of  the  dynamic 
pressure  signals  obtained  on  the  cone  at  Mach  number  6.  At  the  model  base  station  (Sta  40), 
at  the  minimum  free-stream  Reynolds  number  where  transition  was  in  close  proximity  to  the 
model  base,  two  distinguishable  nodes  (like  the  two  humps  on  a  Bactrian  camel)  appeared  in 
the  pressure  fluctuation  spectrum.  Increasing  the  Reynolds  number  from  0.5  x  106  to  1.0  x 
106  produced  an  abrupt  increase  in  the  magnitude  of  the  spectra,  and  the  nodes  reduced 
from  two  to  one  (like  the  hump  on  a  dromedary  camel).  A  further  increase  in  the  free-stream 
Reynolds  number  caused  the  nodes  to  disappear,  and  the  overall  level  of  the  signal 
diminished  slightly.  A  similar  phenomenon  was  observed  in  some  recent  results  obtained  on 
an  8-deg  cone  in  Tunnel  B  at  Mach  number  8  and  is  reported  in  Ref.  37. 

The  variation  in  the  level  of  the  spectra  as  transition  moves  over  the  gage  is  quite  normal; 
that  is,  the  maximum  level  occurs  with  transition  over  the  gage,  the  minimum  value  occurs 
when  the  flow  is  laminar,  and  the  turbulent  condition  produces  a  signal  level  which  is  slightly 
less  than  the  level  existing  when  the  boundary  layer  was  transitional. 

No  unusually  prominent  nodes  appear  in  the  spectra  of  the  small  gage  (0.10-in.  diam)  at 
the  mid-body  station  of  25  in.  at  the  zero  roll  location,  but  the  larger  gage  (0.20-in.  diam)  at 
the  same  station  (25.0  in.)al  the  1 80-deg  roll  location  has  a  single  hump  (peak  at  1 10  kHz)  in 
the  spectra  at  the  minimum  free-stream  Reynolds  number  (0.5  million  per  foot). 

Also,  both  transducers  at  station  25  contain  a  discrete  step  in  the  spectra  at  a  frequency 
of  33  to  35  kHz.  This  type  of  disturbance  was  observed  in  the  surface  pressure  fluctuation 
spectra  of  a  sharp  nose  cone  at  free-stream  Mach  numbers  of  6  and  8  and  is  reported  in  Ref. 

5.  The  source  of  the  disturbances  reported  in  Ref.  5  occurred  at  a  frequency  of  45  to  48  kHz 
and  was  not  identified.  A  similar  disturbance  was  noted  at  56  kHz  in  the  spectra  shown  in 
Fig.  33d  for  pressure  fluctuations  obtained  at  Sta  11.0  in. 

The  results  obtained  in  the  presence  of  a  laminar  boundary  layer  (the  data  from  Fig.  33d)  , 
are  presented  in  terms  of  the  normalized  parameters  (PVq^Uoo/S*  and  wave  number 
fSVUoe,  in  Fig.  34.  This  figure  contains  a  listing  of  the  overall  bandwidth  of  the  fluctuation 
which  includes  95  percent  of  the  energy  of  the  fluctuation,  designated  as  NED  (0.95).  This 
tabulation  simply  shows  that  as  the  free-stream  Reynolds  number  decreases,  the  bandwidth 
of  the  fluctuation  increases. 
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6.4  SUMMARY  OF  TUNNEL  A  AND  B  PRESSURE  FLUCTUATION  RESULTS 

A  summary  of  the  pressure  fluctuation  levels  inferred  or  measured  in  Tunnels  A  and  B  is 
given  in  Fig.  35.  The  band  associated  with  results  inferred  from  the  dynamic  pitot  probe 
reflect  primarily  the  variations  produced  by  the  free-stream  Reynolds  number  variations.  In 
general,  the  FY79  probe  results  are  smaller  than  those  obtained  in  FY78.  The  nozzle  wall 
static  pressure  fluctuations  measured  about  125  in.  upstream  of  the  test  section  are 
consistently  higher  than  the  inferred  free-stream  values.  All  the  results  suggest  that  the 
minimum  pressure  fluctuation  level  occurs  around  Mach  number  4  in  Tunnel  A.  The 
equivalent  sound  pressure  level  (SPL)  in  the  test  section  can  be  computed  as  follows: 

SI’I,  -  201og10C.-  v  20log)0q^  ..  I70.7r,  (52) 

where  q«  is  given  in  lbf/in.2. 

The  curve,  based  on  analysis  from  Ref.  15  and  included  in  the  figure  for  reference 
purposes  only,  is  the  predicted  variation  in  the  pressure  fluctuation  as  a  function  of  free- 
stream  Mach  number  as  formulated  by  Lowson  (sec  Ref.  15).  This  curve  provides  a  rough 
estimate  of  the  Mach  number  trend  in  the  pressure  fluctuation  levels  in  Tunnel  A.  Actually, 
the  hot-wire  results  in  most  cases  (Fig.  19)  suggest  that  the  static  pressure  fluctuations  in 
Tunnel  B  are  less  than  the  values  inferred  from  the  pitot  acoustic  probe  and  also  less  than  the 
values  based  on  the  Lowson  curve. 

The  results  in  Fig.  35b  describe  the  typical  pressure  fluctuation  levels  detected  on  a  body 
of  revolution.  As  expected,  the  free-stream  Reynolds  number  and  the  nature  of  the 
boundary  layer  on  the  body  have  a  direct  effect  on  the  response  of  a  surface  flush-mounted 
dynamic  pressure  transducer.  The  maximum  fluctuations  sensed  by  the  dynamic  pressure 
sensor  on  the  cone  occur  when  transition  is  at  the  gage  location.  Also  as  noted  previously, 
the  gage  size  can  have  a  pronounced  effect  on  the  gage  response:  the  larger  the  gage,  the 
weaker  the  response. 

Figure  36  summarizes  the  attempts  to  identify  the  test  section  static  pressure  fluctuation 
levels  in  Tunnels  A  and  B,  In  general,  in  the  Reynolds  number  range  from  2  to  4  million  per 
fool,  the  static  pressure  fluctuation  levels  deduced  from  the  hot-wire  and  pilot  acoustic 
probes  are  in  reasonably  good  agreement.  The  cone  surface  static  pressure  fluctuations 
increase  with  the  free-stream  Mach  number  and  deviate  significantly  from  the  hot-wire  and 
pitot  probe  results.  Apparently,  there  is  some  amplification  of  (he  free-stream  disturbance  as 
it  propagates  through  the  laminar  boundary-layer  flow  to  the  transducer  on  the  cone 
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surface,  and  this  ainplificaiion  increases  with  the  free-stream  Mach  number.  These  results 
also  suggest  that  the  magnitude  of  the  pressure  fluctuations  decreases  as  the  cone  half-angle 
increases. 

The  correlation  term  used  to  curve  fit  the  free-stream  static  pressure  fluctuation  variation 
with  Mach  number  was  suggested  in  Ref.  15.  The  term  in  the  denominator  is  simply  the 
Eckert  reference  temperature  ratio  with  Tw  being  the  nominal  tunnel  wall  temperature  and 
the  recovery  factor  (77)  equal  to  0.891.  This  curve  fit  was  used  to  estimate  the  pressure 
fluctuation  level  in  Tunnel  C  at  Mach  number  10. 

7.0  CONCLUDING  REMARKS 

The  free-stream  fluctuation  or  turbulence  levels  which  have  been  determined  in  Tunnels 
A  and  B  are  based  on  free-stream  hot-wire  anemometer,  pitot-acoustic  probe,  and  cone 
surface  pressure  fluctuation  measurements.  The  accomplishments  and  significant  results  of 
this  investigation  are  as  follows: 

1.  Reasonable  correlations  based  on  the  analysis  from  Ref.  34  were  established 

# 

between  the  pitot-acoustic  probe  data  and  hot-wire  data.  In  addition  it  was  noted 
that  as  the  free-stream  Mach  number  increased,  the  surface  static  pressure 
fluctuations  measured  on  a  cone  underneath  a  laminar  boundary  layer  deviated 
from  the  free-stream  pressure  fluctuation  levels  based  on  the  hot-wire  and  pitot- 
acoustic  probes.  These  results  should  resolve  most  of  the  discrepancies  that  had 
previously  existed  between  the  AEDC  dynamic  pressure  measurements  and  the 
free-stream  hot-wire  results  noted  previously  in  Refs.  3,  5,  6,  and  7. 

2.  Pressure  fluctuations  measured  on  the  surface  of  a  cone  model  in  the  presence  of 
a  laminar  boundary  layer  may  not  be  a  true  measure  of  the  absolute  level  of  the 
free-stream  pressure  fluctuations  which  characterize  the  free-stream  flow.  This 
fact  should  be  taken  into  account  for  wind  tunnel  tests  in  which  the  influence  of 
free-siream  flow  fluctuations  is  significant,  as  in  the  case  of  naturally  occurring 
boundary-layer  transition. 

3.  In  general,  at  nominal  free-stream  Reynolds  numbers  of  2  to  4  million  per  foot, 
the  free-stream  pressure  fluctuations  as  a  function  of  nozzle  Mach  number, 
tunnel  wall  temperature,  lest  section  static  temperature,  and  dynamic  pressure 
vary  as  follows: 

0.0029  ^/Vf7TJ 
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where 

erVrj  -  i0.5[<Tw/TJ  +  1]  +  0.22tj[(}'  -  1)/2]m2} 

4.  The  results  at  Mach  numbers  2,  4,  and  5  in  Tunnel  A  and  the  pitot -acoustic  data 
in  Tunnel  B  at  Mach  number  6  suggest  that  the  free-stream  fluctuation  ratio  (p/p) 
decreases  with  increasing  free-stream  Reynolds  number;  that  is,  (p/p)  is 
proportional  to  (REde)-°-2S. 

5.  An  evaluation  of  the  pressure  fluctuations  sensed  by  the  pitot-acoustic  probe 
could  not  be  established  until  the  spectral  content  of  the  probe  transducer 
response  signal  was  analyzed.  The  probe  mounting  or  support  system,  the 
transducer  mounting,  and  the  dynamic  calibration  of  the  transducers  could  and  in 
some  cases  did  affect  the  evaluation  of  the  pitot  probe  pressure  fluctuations. 
Basically,  the  spectral  distribution  provides  the  method  for  identifying  those 
dynamic  disturbances  sensed  by  the  probe  which  are  not  a  part  of  the  test  section 
sound  source  (turbulence). 

6.  A  by-product  of  this  study  was  the  development  of  the  software  and  analysis  used 
to  interface  the  SD360  signal  processor  with  the  facility  computer.  The  signal 
processor  was  used  to  convert  the  analog  signals  into  a  digital  spectral  format 
(transform,  power  spectral,  or  cross-spectral  density),  and  the  facility  computer 
was  used  to  normalize  and  plot  the  spectra,  digitally  filter  and  recompute  the 
root -mean-square  of  the  signals,  and  provide  additional  digital  signal  processing 
for  further  analysis  of  the  signal  results. 

7.  The  stilling  chamber  pressure  fluctuations  were  found  to  correlate  simply  as  a 
function  of  the  mean  value  of  the  stilling  chamber  pressure.  As  expected,  the 
addition  of  high-pressure  makeup  air  to  maintain  the  flow  rate  in  the  tunnel 
circuit  was  found  to  significantly  increase  the  stilling  chamber  turbulence.  The 
hot-wire  measurements  obtained  in  1965  and  the  present  surface  static  pressure 
fluctuation  measurements  were  in  good  agreement. 

In  addition  to  these  accomplishments,  this  present  effort  revealed  the  following 

problems:  - 

1.  Additional  test  results  and  analysis  are  needed  at  Mach  numbers  6  and  8  to 
understand  the  trend  in  which  the  pressure  fluctuation  ratio  based  on  the  hot-wire 
measurements  increases  with  the  free-stream  Reynolds  number. 
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2.  In  the  present  studies,  the  fluctuation  measurements  were  made  at  one  point  in 
the  test  section.  The  question  remains  concerning  the  spatial  variation  in  the 
fluctuations  throughout  the  test  section.  In  particular,  the  presence  of  a  small 
region  of  nonuni  form  flow  along  the  axis  of  the  axisymmetric  nozzles  of  Tunnel  B 
suggests  that  some  effort  is  warranted  to  define  the  radial  variation  in  the  tunnel 
test  section  turbulence.  Since  the  sound  disturbances  emanate  from  the  tunnel 
wall  turbulent  boundary  layer,  there  is  a  need  to  map  the  variation  in  the 
turbulence  level  near  the  tunnel  wall. 

3.  The  possibility  exists  for  gleaning  additional  information  from  the  present  flow 
quality  turbulence  data  by  examining  the  correlations  that  may  exist  between  the 
various  types  of  measurements.  This  additional  analysis  would  require  some 
additional  programming  and  evaluation  of  the  existing  data  source.  Possibly,  the 
convective  velocity  of  the  sound  source  associated  with  the  flow  over  a  typical 
conical  body  could  be  identified  for  various  wind  tunnel  conditions. 

4.  The  turbulence  levels  sensed  by  the  instrumentation  in  a  planar  surface  aligned 
with  the  flow  need  to  be  established.  This  would  require  reinstrumenting  the  flat 
plate  and  testing  it  in  Tunnels  A  and  B. 

5.  The  results  presented  in  this  report  should  be  supplemented  by  an  investigation  of 
the  flow  fluctuations  which  characterize  the  interior  of  the  boundary  layer  above 
the  surface-mounted  transducers.  Such  an  effort  would  address  the  need  for  an 
improved  understanding  of  the  differences  between  free-stream  and  surface 
pressure  fluctuations.  A  limited  inquiry  for  this  purpose  could  be  made  using 
existing  hot-wire  anemometer  measurements  acquired  during  the  present  test  and 
earlier  tests  using  slender  cone  models  in  Tunnel  B.  However,  the  reduction  and 
analysis  of  these  data  for  this  purpose  were  beyond  the  scope  and  resources  of  the 
earlier  tests  and  of  the  present  effort. 
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Fig.  1  Tunnel  A 
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c.  Upstream  of  Tunnel  A  nozzle 
Figure  1.  Continued. 
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d.  Stilling  Chamber  (as  of  FY78) 
Figure  1.  Continued. 
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e.  Transducer  installation  in  the  stilling  chamber 
Figure  1.  Continued. 
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g.  Nozzle  —  diffuser  section 
Figure  1.  Concluded. 


AEDC-TR-01  -5 


Figure  2.  Typical  installation  of  the  tunnel  wall  flush-mounted  dynamic  pressure  sensors. 
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Figure  3.  Tunnel  A  flat  plate  installation  (FY78). 
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Figure  4.  Tunnel  A  installation  with  nominal  locations  of  the  dynamic  pressure  sensors  (transducers). 
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Figure  6.  Test  installation  in  Tunnel  A  (FY79). 
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b.  Tunnel  test  section 
Figure  7.  Tunnel  B. 
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a.  Installation  sketch 
Figure  8.  Test  installation  in  Tunnel  B. 
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Figure  10.  Flat  plate  installation  in  Tunnel  B,  typical  shadowgram  at  Mach  number  8 
and  Reynolds  number  2  x  106  per  feet. 
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1 


68 


AEDC-TH  -81  -5 


Exposed  Transducer 


Vented  to 
Model  Cavity 


Blind  Transducer  Model 
Wall  Installation 


b.  Installation  codes 
Figure  11.  Concluded. 
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a.  Details 

Figure  12.  Hot-wire  anemometer  probe. 
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b.  Silhouette  photograph 
Figure  12.  Concluded. 
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Figure  15.  Tunnel  A  stilling  chamber  pressure  fluctuations. 


NOTE:  Open  Symbols  ^  Nov.  1965  Hot  Wire  Data  Analysis 

Solid  Symbols  ^  FY  1978  Wall  Static  Pressure 
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b.  Velocity  flunctuations 

Figure  16.  Summary  of  stilling  chamber  turbulence  measurements  in  Tunnel  A. 
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c.  PT  =  27.0  psia,  PT/PT  =  0.011  percent,  Af  =  100  Hz,  three-stage  plant  configuration. 

Figure  17.  Continued. 


d.  PT  91.3  psia,  PT/PT  0.0020  percent,  Af  -  100  Hz,  four~stage  plant  configuration. 

Figure  17.  Continued. 
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f.  FT  =  137.2  psia  (high-pressure  makeup  air  added),  PT/PT  =  0.0462  percent, 
Af  =  200,  four-stage  plant  configuration. 

Figure  17.  Concluded. 
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Pressure  Fluctuation,  P(f)/P  (f  =  150  Hz) 
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h.  Normalized  velocity  fluctuations 
Figure  22.  Concluded. 
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Pressure  Fluctuation,  P{t)/P  (f  =  600  Hz) 
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b.  Normalized  velocity  fluctuations 
Figure  23.  Concluded. 
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a.  Normalized  pressure  fluctuations 

Figure  24.  Normalized  distributions  of  fluctuations  inferred  from  hot-wire  measurements, 
M,  =  8,  RE.  =  3  x  106  per  foot. 
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b.  Normalized  velocity  fluctuations 
Figure  24.  Continued. 
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Temperature  Fluctuation,  T(i)/T  (f  =  600  Hz) 
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=  4.0,  RE  =  1.0  x  106  per  foot 
Figure  26.  Concluded. 
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Figure  28.  Tunnel  A  tare  for  the  dynamic  pressure  sensor  of  the  pitot  acoustic  probe  (FY78). 
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Figure  29.  Tunnel  B  tare  for  the  dynamic  pressure  sensor  of  the  pitot  acoustic  probe  (FY78). 
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a.  Normalized  energy  distribution 

Figure  30.  Typical  spectral  distribution  in  the  wall  static  pressure  at  tunnel  Station  125  in. 
(Mach  number  4.0,  RE.,  =  4  x  106  per  foot). 
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Figure  32.  Continued. 
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a.  PCB®113  at  Station  40.0  (3.33  ft) 

Figure  33.  Pressure  fluctuation  spectra  on  a  6-deg  sharp  nose  cone  at 
Mach  number  6.0  (zero  angle  of  attack). 
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c.  PCB®113  at  Station  25.0  (2.08  ft},  0  =  180  deg 
Figure  33.  Continued. 
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d.  PC  B®  1.5  at  Station  11.0  (0.92  ft) 
Figure  33.  Concluded. 
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Figure  34.  Normalized  power  spectral  density  of  the  surface  pressure  fluctuations  on  a  6-deg  sharp  nose  cone. 
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a.  Inferred  static  pressure  fluctuations  based 
on  dynamic  pitot  pressure  probe 


b.  Typical  cone  (6  deg)  surface  static  pressure  fluctuation 
levels  obtained  with  a  O.IO-in.-diam  PCB® 

Figure  35.  Typical  pressure  fluctuation  levels  in  the  AEDC  Wind  Tunnels  A  and  B. 
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Figure  36-  Summary  of  pressure  fluctuation  levels  in  AEDC  Tunnels  A,  B  and  C 
at  nominal  RE„  =  2  x  106  to  4  x  106  per  foot. 
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APPENDIX  A 

HOT-WIRE  ANEMOMETER  CALIBRATION 

Two  types  of  calibration  are  required  for  the  hot-wire  probes  if  quantitative  data  are  to 
be  extracted  from  flow-field  measurements  made  using  hot-wire  anemometry  techniques. 


The  first  calibration  is  made  to  determine  certain  physical  characteristics  of  the  sensor: 
wire  resistance  at  zero  current  at  a  reference  temperature  Tf  =  32° F,  the  first-order 
temperature  resistivity  coefficients  at  the  reference  temperature,  and  the  wire  effective 
length.  This  calibration  is  done  in  the  VKF  instrumentation  laboratory  prior  to  the  wind 
tunnel  measurements,  and  a  posttest  calibration  is  performed  on  those  wires  which  survive 
wind  tunnel  use.  Probes  are  placed  in  a  laboratory  heat-treating  furnace,  and  the  sensor 
resistance  is  measured  at  a  series  of  six  temperature  levels  from  room  temperature  up  to 
500°F  for  wires  used  in  Tunnel  A  and  seven  temperature  levels  up  to  900°F  for  wires  used  in 
Tunnel  B.  At  each  oven  temperature  level  the  sensor  resistance  is  measured  at  each  of  six 
wire  heating  currents  up  to  1 .2  ma  for  50-jdn.-diam  wires  and  up  to  0.6  ma  for  20-/nn.-diam 
wires.  At  each  temperature  level  the  results,  in  the  form  of  measured  resistance  versus 
heating  current-squared  (for  example,  see  Fig.  A-l),  are  extrapolated  to  zero  current  using  a 
first  order,  least  squares  curve  fit.  Data  from  the  various  temperature  levels  are  combined  in 
the  form  of  zero-current  resistance  versus  temperature  (Fig.  A-2)  according  to  the  formula 
R  =  Rr[l  +  «r(T  -  Tr)]^  A  first-order,  least  squares  fit  of  the  data  yields  the  sensor 
reference  resistance,  Rr,  at  T  =  Tr  and  the  linear  coefficient  of  resistivity.  Wires  of  platinum 
— 10  percent  rhodium  alloy  of  10-jtin.  diam  were  found  to  have  a  zero  second-order 
temperature  coefficient  of  resistivity  by  Philco-Ford  (Ref.  22)  and  earlier,  for  SO-fiin.-diam 
wires,  by  Spangenberg  (Ref.  39).  This  finding  was  assumed  to  be  valid  for  the  20-  and 
50-jiin.-diam  wires  used  for  the  tunnel  free-stream  measurements  of  the  present  report. 

p 

The  second  calibration  of  the  probes  is  made  in  the  wind  tunnel  to  determine  the 
relationship  of  wire  heat-loss  coefficient  (Nusselt  number)  and  wire  temperature  recovery 
factor  as  functions  of  the  wire  Reynolds  number.  (For  example,  see  Fig.  A-3.)  This  Reynolds 
number  (Ref.  29)  is  calculated  using  local  density  and  velocity  and  the  dynamic  viscosity 
based  on  local  total  temperature.  The  wire  diameter  is  used  as  the  characteristic  length.  The 
Reynolds  number  variations  are  obtained  by  varying  the  tunnel  total  pressure,  PT,  while 
holding  the  tunnel  total  temperature,  TT,  nominally  constant.  Each  wire  was  calibrated  at 
the  free-stream  Mach  number(s)  for  which  it  was  used.  The  Nusselt  number  and  recovery 
factor  can  be  considered  independent  of  Mach  number  for  M«,  >  1.5  (Ref.  29). 
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Squared  Wire  Heat 

Figure  A-1.  Typical 
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Figure  A-3.  Typical  flow  calibration  results. 
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APPENDIX  H 

DYNAMIC  PRESSURE  SENSOR  CALIBRATION 
B.1.0  INTRODUCTION 

In  some  cases,  dynamic  calibrations  of  the  DPS  transducers  were  obtained  using  the  VKF 
shock  tube.  In  each  calibration,  the  shock  tube  subjected  the  transducer  to  a  known 
transient  pressure  input  closely  resembling  a  step.  Analysis  of  the  transducer  response 
generated  a  useful  calibration  in  the  form  of  a  transfer  function.  A  description  of  the 
apparatus  and  general  procedure,  with  results,  follows. 

Tables  B-l  and  B-2  summarize  the  d-c  calibration  results  of  the  transducers  used  in  this 
study.  In  some  cases  the  effects  of  subatmospheric  environmental  pressure  level  on  the  gage 
calibrations  are  included. 


B.2.0  APPARATUS 


B.2.1  SHOCK  TUBE 

The  VKF  shock  tube,  shown  in  Fig.  B-l,  is  approximately  10  fi  long  with  an  inside 
diameter  of  2,5  in.  A  Mylar®  diaphragm  separates  the  high-pressure  “driving”  section  from 
the  low-pressure  “driven"  section,  the  latter  being  approximately  8  ft  in  length.  When  a 
sutficiem  pressure  differential  is  developed  across  the  diaphragm,  it  ruptures,  producing  a 
shock  wave.  The  wave  travels  down  the  tube  past  a  strut-supported  hemisphere-cylinder 
pressure  probe.  The  transducer,  mounted  flush  at  the  probe  tip,  experiences  a  sudden  rise  in 
pressure  and  temperature  corresponding  to  the  establishment  of  supersonic  flow  behind  the 
shock.  Following  400  to  850  ^sec  of  steady  flow  the  reflected  waves  return  and  result  in  flow- 
breakdown. 

The  shock  wave  strength  and  velocity  determine  the  magnitude  of  the  pressure  rise 
imposed  on  the  probe.  These  factors  are,  in  lurn,  determined  by  the  pressure  ratio  across  the 
diaphragm.  As  the  shock  tube  lacks  a  mechanical  diaphragm  puncture  mechanism,  the 
primary  controls  arc  the  driven  section  charge  pressure  and  diaphragm  thickness.  To  fire  the 
tube,  the  driver  section  pressure  is  gradually  increased  until  the  diaphragm  ruptures. 
However,  as  will  be  shown  in  Section  B.4.0,  the  shock  tube  conditions  are  highly  repeatable. 
The  measured  parameters  include  the  driven  section  temperature,  the  driven  section  initial 
pressure,  and  Ihe  shock  velocity.  The  shock  velocity,  (Vs),  is  determined  by  measuring  the 
time  required  for  the  shock  to  travel  a  known  distance  along  the  tube.  Two  highly  responsive 
thin  film  resistance  thermometers,  located  1.001  ft  apart  along  the  lube,  are  used  to  detect 
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the  passing  of  the  shock.  An  electronic  counter  measures  the  time  between  the  responses  of 
the  temperature  probes  (see  Fig.  B-2). 

B.2.2  DATA  SYSTEM 

The  data  acquisition  hardware  is  illustrated  schematically  in  Fig.  B-l.  During  a 
calibration  run  the  transducer  output  signal  is  conditioned,  amplified,  and  recorded  on 
analog  tape  at  maximum  tape  speed,  120  in. /sec.  Following  the  run,  the  spectral  analyzer  is 
used  to  digitize  the  analog  data.  To  maximize  the  number  of  digital  data  points 
characterizing  the  short  event,  the  tape  is  run  at  a  lower  speed.  30  in. /sec,  during  playback 
while  the  analyzer  samples  at  a  maximum  frequency  of  307.2  kHz.  Compensating  for  the 
factor  of  four  reduction  in  tape  speed  produces  an  effective  sample  frequency  of  1228.8 
kHz.  The  analyzer  utilizes  1,024  points  to  define  the  transducer  output.  It  is  capable  of 
fiUering  the  output  to  prevent  aliasing,  transforming  the  output  into  the  frequency  domain, 
and  obtaining  transfer  functions  with  the  assumption  of  a  pure  step  input.  The 
minicomputer  records  the  analyzer  results  on  disk.  The  data  are  then  available  for  plotting 
on  a  graphics  terminal.  The  minicomputer  also  serves  as  an  interface  with  the  facility 
computer.  The  calibrations  can  be  transferred  to  the  facility  computer  for  use  in  more 
detailed  calibration  analysis  and  in  wind  tunnel  test  data  reduction. 

B.3.0  DATA  REDUCTION 

A  transfer  function  is  generated  by  analyzing  the  transducer  response  to  a  known 
transient  pressure  input.  The  pressure  input  produced  in  the  shock  tube  can  be  characterized 
as  a  perfect  step  or  a  modified  step,  technically  representing  all  frequencies  in  the  spectrum. 
The  modified  step,  the  overshoot,  accounts  for  the  transient  establishment  of  the  probe  bow 
shock,  requiring  a  finite  time  following  the  passage  of  the  tube  shock.  Modified  inputs  can 
be  generated  analytically  or  by  using  extremely  responsive  transducers,  such  as  the  PCB-1 13, 
to  define  the  shape  for  less  responsive  transducers.  Pure  step  inputs  are  utilized  in  this 
Appendix  for  illustration  purposes.  Figure  B-2  presents  the  analytical  expressions  used  to 
calculate  the  step  input  from  shock  tube  measurements. 

The  analysis  of  the  dynamic  pressure  data  is  performed  in  the  frequency  domain.  The 
time  domain  and  frequency  domain  representations  of  analog  data  are  related  by  the 
continuous  Fourier  transform  pair: 

11(f)  -  f  °°  h(i)e-i2?rfldi  (Forward  Transform) 

•*—oo 

Mi)  -  P  iKOeJ^r,  df  (Inverse  Transform) 
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Digital  (sampled)  data  are  transformed  using  the  discrete  Fourier  transform: 

N— i 

H/-!LA  =  £  Mkt)e  (Forward  Transform) 

^  k=o 


h(kt)  =  — 
N 


eJ27rnk/N 


(Inverse  Transform) 


The  spectral  analyzer  and  facility  computer  calculate  the  discrete  transform  using  an 
algorithm  called  the  fast  Fourier  transform,  or  FFT.  The  FFT  greatly  reduces  the  number  of 
numerical  operations  required  to  define  the  transform,  and  therefore  the  FFT  is  a  practical 
tool  for  dynamic  data  analysis. 


A  data  point  in  the  frequency  domain  can  be  represented  by  a  complex  number  or  an 
amplitude  and  phase  angle.  The  representations  are  related  as  follows; 

|H(0|  =  +  I2(f) 


0.(0  =  Tan'HKfVmD] 

where 


|H(f);  =  Amplitude  of  the  transform 
(1>( f)  =  Phase  angle 

R(f)  =  Real  part  of  the  complex  representation 
1(1)  =  Imaginary  part 

The  amplitude/phase  angle  representation  was  used  as  the  dynamic  pressure  transducer 
calibration. 

Filtering  is  required  when  the  data  are  digitized  to  prevent  frequency  foldover,  or 
aliasing,  in  the  frequency  domain.  The  sample  rate  must  be  at  least  twice  the  maximum 
frequency  present  in  the  analog  signal.  The  spectral  analyzer  applies  the  appropriate  low- 
pass  filter  when  the  data  are  digitized.  The  filter  cutoff  frequency  of  the  analyzer, 
corresponding  to  the  307.2-kHz  maximum  sample  rate,  is  120  kHz.  When  one  compensates 
for  the  tape  speed  reduction  (a  factor  of  4)  during  digitization,  the  effective  cutoff  frequency 
becomes  480  kHz. 
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Leakage  in  the  frequency  domain  is  reduced  by  applying  one  of  several  types  of 
windows.  The  spectral  analyzer  windows  the  data  using  either  a  weighted  rectangular 
window  or  a  Kaiser-Bessel  weighting  function.  The  calibration  program,  on  the  facility 
computer,  uses  a  subtraction  scheme  to  window  the  data.  A  point  near  the  center  of  the 
data,  point  570  for  example,  is  selected  as  the  start  point  of  the  window.  The  remaining 
points  are  then  redefined  by  subtracting  consecutive  points  in  the  first  half  of  the  data  from 
the  window  start  point.  The  result  is  that  the  data  contain  the  transient  characteristics  of  the 
original  signal  with  equal-valued  end  points.  Other  window  techniques,  such  as  the  Kaiser- 
Bessel  weighting,  can  easily  be  added  to  the  calibration  program.  The  effects  of  the  window 
cancel  during  the  calculation  of  the  transfer  function. 

The  procedure  for  generating  pressure  transducer  transfer  functions  on  the  facility 
computer  is  illustrated  in  Fig.  B-3.  The  transducer  output  signal  and  the  input  curve  are 
compared  to  determine  the  d-c  scale  factor.  The  transducer  output  signal  is  multiplied  by  the 
d-c  scale  factor  to  convert  from  mv  to  psi.  The  input  curve  and  output  signal,  in  psi,  are 
windowed  and  plotted  together  to  check  for  time  alignment.  A  misalignment  is  corrected  by 
time  shifting  (signal  rotation).  The  input  curve  and  output  signal  are  then  transformed  into 
the  frequency  domain.  The  transfer  function  is  formed  by  calculating  the  quotient  of  the 
amplitudes  and  the  difference  between  the  phase  angles  as  follows: 

Transfer  ^^InpUl  “  ^^Oulput 

The  transfer  function,  applied  to  wind  tunnel  lest  data,  removes  transducer  response 
effects  and  thus  reveals  the  true  pressure  input.  The  procedure  is  presented  in  Fig.  B-4.  The 
transducer  output  is  converted  to  psi  through  multiplication  with  the  d-c  scale  factor  and  is 
transformed  into  the  frequency  domain.  The  transfer  function  is  applied  to  obtain  the  true 
input  to  the  transducer: 

lHU)lT..t  -|H<0,T,„.(.r-|H«,  Teat 

Inp'Jt  Output 

Test"  ^ Transfer  +  Test 
Input  Output 


The  true  input,  representing  the  pressure  fluctuations  experienced  by  the  transducer,  is  then 
ready  for  analysis. 
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B.4.0  RESULTS  AND  DISCUSSION 


The  following  dynamic  pressure  transducers  were  calibrated  in  support  of  the  present 
research  effort: 


Transducer 

Serial  Ho 

PCB-1 13A21 

1446 

PCB-1 13A21 

1447 

PCB-1 12A21 

1407 

PCB-1 05A 

568 

PCB-1 06A 

248 

Kulite  CQ-080-5 

J1-36 

Kulite  CQ-080-5 

J1-37 

Endevco  8510-5 

DM- 43 

Endevco  8510-5 

DN-76 

Three  optimum  shot  conditions  were  used: 


Diaphragm  T1 , 

,  °R 

PI,  psi 

4t,  usee 

PD1-P1,  psi 

2  Layers,  0.002-in.  Mylar  533 

0.327 

321 

7.40 

1  Layer,  0.002-in.  Mylar 

0.135 

312 

3.  36 

1  Layer,  0.002-in.  Mylar 

0.0234 

262 

1.02 

The  repeatability  in  the  shot  conditions  is  illustrated  in  Fig.  B-5. 

Sample  calibrations  are  presented  in  Figs.  B-6  through  B-9  for  the  PCB-1 12  A2I  (1407), 
PCB-1 13A21  (1447),  PCB-105A  (56S),  and  the  ENDEVCO  B510-5  (DN-43),  respectively. 
Each  Figure  presents  the  example,  step  by  step,  in  accordance  with  Fig.  B-3.  It  should  be 
noted  that  the  signal  noise  was  not  averaged  out  in  the  examples. 

The  PCB-I12A21  and  PCB-105A  have  manufacturer’s  quoted  resonant  frequencies  of 
250  kHz  which  are  clearly  visible  in  the  transfer  functions.  The  response  is  relatively  flat  at 
the  frequencies  below  resonance.  The  PCB-1 13A21  has  a  manufacturer’s  quoted  resonant 
frequency  of  500  kHz,  which  is  beyond  the  useful  frequency  range  of  the  transfer  function. 
The  decay  in  the  PCB-105A  transducer  output  is  due  to  the  temperature  sensitivity  of  the 
transducer,  which  affects  only  low-frequency  response  of  the  transducer.  This  temperature 
effect  has  a  negligible  influence  on  the  overall  dynamic  characteristics  of  the  transducer.  An 
RTV  cap,  nominally  0.015  in.  thick,  placed  over  the  transducer  attenuates  the  overall  level 
of  the  response  about  6  db.  Although  not  shown,  the  variation  of  the  transfer  function 
across  the  frequency  spectrum  is  not  significantly  affected. 
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The  remaining  factor  which  can  alter  the  output  of  the  transducer  is  its  physical  size  and 
the  local  convective  velocity  of  the  pressure  disturbance  over  the  transducer  surface.  The 
diameters  of  the  present  transducer  varied  in  size  from  0.099  to  0.40  in.  Most  of  the  results 
presented  in  this  report  were  obtained  with  a  0.099-in. -diam  transducer. 

G.  M.  Corcos  (Ref. 35)  formulated  a  theoretical  correction  factor  to  account  for  the 
attenuation  in  the  output  of  a  finite  size  transducer.  For  convience  purposes  and  future 
reference,  these  correction  factors  were  curve  fitted  as  follows. 


1.  Range:  0.050  >  ojr/ut  >  1.00 


Ao  =  3.01394 

A,  =  -4.33121 
A2  =  3.63664 

A3  =  -1.68194 
A4  =  0.39941 

A;  =  -0.03814 

2.  Range:  1.00  >  cur/uc  >  2.40 


A0  - 

1.25246 

A|  = 

-0.48254 

A2  = 

0.098501 

a3  = 

-0.011444 

a4  = 

0.00070593 

As  = 

-0.000017893 

3.  Range:  2.40  >  cdr/uc  a  3.5 

f  /F  =  3.29364  exp  1-1.50675 (tur/u  )] 

m  c 


4.  Range:  3.5  >  wr/uL. 


F'F 


B  (<5>r/u  )  Bl 
0  r 
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Range  BQ  B1 

» 

3.5  >  a) r/u  >  5.0  0.090844  1.3184 

c  — 

5.0  >  ur/u  0.50127  2.3835 

c 

Figure  B-10  illustrates  the  type  of  correction  that  would  have  to  be  applied  to,  for 
example,  a  power  spectrum  for  each  of  the  three  sizes  of  transducers  used  in  the  present 
study  assuming  the  local  convective  velocity  is  1,140  ft/sec.  Since  a  major  portion  of  the 
present  data  was  obtained  with  a  0.099  in.-diam  PCB,  a  smaller  transducer,  and  in  a  flow 
field  with  higher  convective  velocities;  the  Corcos  attenuation  factor  was  not  applied  to  the 
present  results.  Note  that  the  Corcos  factor  is  at  a  maximum  at  the  higher  frequencies  where 
the  signal  content  in  the  measured  pressure  spectra  is  usually  smallest  (i.e.,  negligible). 
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Figure  B-1. 
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To  Tape 
Machine 


ToTape 

Machine 


b.  PCB®  transducer  signal  conditioning 
Figure  B-1.  Continued. 


To  Tape 
Machine 


To  Tape 
Machine 


Transducer 

Excitation,  VOC 

Scale  Factor,  mv/psi 

Bridge  Resistance;  Q 

Jl-36  (Kulite) 

5 

15.7 

2,500 

Jl-37  (Kulite) 

5 

15.5 

2,500 

Gl-51  (Kulite) 

15 

12. 0 

3,500 

8510-5  (Endevco) 

10 

=  70.0 

1,800 

c.  Kulite  and  Endevco  transducer  signal  conditioning 
Figure  B-1.  Concluded. 
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J _ I 

fw—  Al.  -•“]  T’  me 


Calculation  of  Step  Input 

I.  Tube  Shock  Wave  Velocity  and  Mach  Number: 
Us  =  1.001/At 


M  = 
s 


=  us/[(49.0223)^] 

II.  Mach  Number  in  Supersonic  Flow  Region  Behind  Shock: 


M2  "  2Vis  ~  1 


2-vM  -  (y 
s 


1)  (v  -  +  2 


III.  Pressure  Ratios  Based  on  Mn  and  M  : 

2  s 


p2/p 


i  -K  ■ 

"p/P2  ‘  [ 


(V  -  1) 


/ 


(Y  +  1) 


-i  Y/(Y-D 


(  y  +  i  m2/2 


x  Jty  +  1)  j( 2 


V  +  1)  /(2yM“  -  (Y 


”)] 


1/(Y-1) 


TV.  Pressure  Srep. 

Pp  -  (VpJ(vri)(Pi> 

Figure  B-2.  Calculation  of  shock  tube  pressure  step. 
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A llg ned/ Windowed  Pressure 


Transform  of  Input 


Freq,  kHz 


Figure  B-3.  Dynamic  calibration  data  reduction. 
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Time,  psec 


NOTE:  Phase  Angle  Plots 

Not  Shown  I  Multiply  by  DC 

?  Scale  Factor 


Output, 

psf 


Transducer  Output,  psl 


Time,  psec 


Transfer  Function 
from  Calibration 


Mag 

I  nput  /Output 


Frequency 


I 

Transform  of  Output 


Multiply:  W  x  Output  ■  I  nput 


Figure  B-4.  Wind  tunnel  test  data  reduction  using  dynamic  calibration. 
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Input  Transform,  psIJHz 
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d.  Transform  of  step  input 
Figure  B-6.  Continued. 
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Output  Transform,  pst/Hz 
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f.  Transfer  function 
Figure  B-6.  Concluded. 
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Transducer  Output,  mv 
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0  80  160  240  320  400  480  560  MO  720  800  880 

Time,  usee 

a.  Step  input 


Time,  psec 

b.  Transducer  response  to  step  input 
Figure  B-8.  Dynamic  calibration  of  a  PCB-105A. 
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Time,  psec 

c.  Input  and  output  aligned  in  time  and  windowed 

800  t 

720- 

640 

560 

480- 

400- 


d.  Transform  of  step  input 
Figure  B-8.  Continued. 
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a.  Step  input 


-50 - . - - - , - . - , - - - - - , - , 

0  80  160  240  320  400  480  560  640  720  800  880 

'  Tine,  nsec 

b.  Transducer  response  to  step  input 
Figure  B-9.  Dynamic  calibration  of  an  Endevco  8510-5  transducer. 
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720- 


d.  Transform  of  step  input 
Figure  B-9.  Continued. 
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Output  Transform.  pst/Hz 


AEDC-Tfl-ei-S 


Freq,  kHz 


e.  Transform  of  transducer  response 


f.  Transfer  function 
Figure  B-9.  Concluded. 
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Table  B-1.  FY79  Transducer  Calibrations 


Serial 

No. 


Scale  Factors,  mv/psi 


Subatmospheric  /^tory  , 
Calibration  (Atmospheric) 
Calibration 


Location 


Kul ite 

Jl-36 

_ 

f 16 . 80 1  * 
1-2.17/ 

15.7 

Stilling 

Chamber 

Jl-37 

/ 16.391* 
1-2.06) 

15.5 

Nonoperating 

Sidewall 

Gl-51 

/ 13.141* 

1 -0.57) 

12.0 

Operating 

Sidewall 

PCB 

113 

1446 

f  26. 381* 
l -0.33) 

22.0 

Flat  Plate, 
PW1 

PCB 

106 

248 

84.8  i  4.4 

100.0 

Flat  Plate, 
PW2 

PCB 

112 

1407 

)  52.721* 
l  -0 . 78) 

42.0 

Window  Wall , 
PWW 

PCB 

]  05 

580 

28.0  ±  2.0 

30.0 

Tank,  PMIT 

PCB 

105 

568 

25.4 

26.0 

Pitot  Probe 

B&.K 

Hie 

803.2 

N/A 

♦Scale  factor  sensitive  to  local  pressure:  SF  =  mP  + 

Not  enough  data  to  fully  evaluate  Gages  5,  7  through  9  for  scale 
factor-pressure  effects. 
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Table  B-2.  FY80  Transducer  Calibrations 


Gage 

No.» 

Serial 

No. 

Scale  Factors,  mv/psi 

Type 

Subatmospheric 

Calibration 

Factory 

(Atmospheric) 

Calibration 

Location 

1 

PCB 

105 

845 

24.5 

Sta  U.O 

2 

828 

17.2 

Sta  12.0, 
Blind 

3 

846 

28.5 

Sta  11.0 

4 

PCB 

113 

1447 

(  28.04  1** 

\  -0.36  1 

23.0 

Sta  12.5 

5 

PCB 

105 

826 

20.6 

Sta  25.0 

6 

ENDEVCO 

DM-43 

f  68.86  !♦* 

\  -3.64  ) 

N/A 

Sta  25.0, 
Blind 

7 

PCB 

113 

1446 

(  26.38  \** 

\  -0.33  / 

22,0 

Sta  25.0 

8 

ENDEVCO 

DK-76 

I  71.94  )♦♦ 

\  -2.06  J 

N/A 

Sta  25.0 

9 

PCB 

112 

1407 

(  52.72  )♦♦ 
l  -0.78  J 

42.0 

Sta  40.0 

10 

PCB 

105 

829 

18.5 

Active,  Pitot 

1 1 

PCB 

105 

827 

18.2 

Blind,  Pitot 

♦Gages  1  through  9  on  6-deg  cone:  Sta  0  at  cone  nose. 
♦♦Scale  factor  sensitive  to  local  pressure;  SF  =  mP  +  B; 
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The  basic  parameters  characterizing  the  general  features  of  the  tunnel  wall  turbulent 
boundary,  which  contributes  significantly  to  the  test  section  turbulence,  are  Re*,  6*.  and  5. 
These  boundary-layer  parameters  were  extracted  from  Refs.  40  and  41  and  are  summarized 
in  Figs.  C-l  through  C-4.  It  is  interesting  to  note  that  the  Reynolds  number  based  on  the 
tunnel  wall  boundary-layer  momentum  thickness  appears  to  be  solely  a  function  of  the  lest 
section  Reynolds  number  in  Tunnel  A,  which  is  an  uncooled,  two-dimensional  nozzle.  In 
Tunnel  B,  which  is  a  water-cooled  ,  axisymmetric  nozzle,  Re*  is  dependent  upon  the  tunnel 
unit  Reynolds  number  and  nozzle  Mach  number. 


The  Van  Driest  II  analysis  (Ref.  42)  for  thin,  turbulent  boundary  layers  in  compressible 
streams  with  zero  pressure  gradients  was  used  to  estimate  the  growth  of  the  boundary  layer 
along  the  Tunnel  A  nozzle  wall.  In  this  analysis,  the  local  Reynolds  number,  RE„,  was 
evaluated  as  an  integral  function  of  the  local  nozzle  stream  properties;  that  is, 

“W.  (pu'>)x  d\ 


The  local  boundary-layer  profiles  and  thickness  parameters  were  evaluated  at  equal  intervals 
along  the  nozzle  wall  from  the  throat  to  the  test  section  of  Tunnel  A.  A  comparison  of  the 
experimental  values  of  RE*  with  the  theoretical  estimates  is  given  in  Fig.  C-4.  The  theoretical 
results  show  a  Mach  number  effect  which  is  not  present  in  the  experimental  results.  This 
discrepancy  may  be  due  to  the  fact  that  in  the  theoretical  analysis  the  pressure  gradient 
effects  on  the  nozzle  wall  boundary  layer  were  neglected,  the  boundary  layer  was  assumed  to 
be  fully  turbulent  from  the  nozzle  throat  to  the  test  section,  and  the  boundary  layers  were 
quite  thick. 


The  experimentally  defined  values  of  RE*  based  on  results  from  Ref.  40  and  the  Winter 
and  Gaudet  correlation  (Ref.  43)  were  used  to  estimate  the  nominal  test  section  wall  skin 
friction  coefficient  given  in  Fig.  C-5.  In  this  correlation,  which  was  found  to  agree  with  wind 
tunnel  results  presented  in  Ref.  44,  the  skin  friction  is  defined  as  follows: 

ef  “  [O.OltO(FgHE0)-°‘l77S]/'Fr 


where 

and 


Fg  =  i  i  0.056M^ 
Fc  =  (l  i  0.2M2)°-3 
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As  noted  in  Section  2.0  the  local  friction  or  shear  coefficient  can  be  shown  to  be  roughly 
proportional  to  the  test  section  free-stream  pressure  fluctuations.  Therefore,  the  parameter 
(cr)  would  be  a  useful  correlating  (or  normalizing)  parameter  for  free-stream  fluctuation 
measurements. 

Figure  C-6  defines  the  shape  factors  associated  with  the  test  section,  tunnel  wall,  and 
turbulent  boundary-layer  thickness  parameters.  There  is  reasonable  agreement  between  the 
Tunnel  A  boundary-layer  shape  factor  ratio  based  on  “6*/0”  and  the  semi-empirical 
correlation  given  in  Ref.  20  for  other  wind  tunnel  shape  factors. 
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a.  Tunnel  A  side  wall 


wall  momentum  thickness. 
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Tunnel  C: 


TO 


a.  Contoured  nozzle/flexible  plate 
Figure  C-2.  Tunnel  wall  boundary-layer  displacement  thickness. 
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HEm  x  ]0"6  (ft-1) 


b.  Tunnel  A  side  wall 
Figure  C-2.  Concluded. 
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20. 


10. 


c  ,  in . 


10.0 


i,  in. 


1.0»- 
1  .0 


1  .o 


RE,;  X  10-6  .  rt- 

a.  Tunnel  A  side  wall 


LO.O 


b.  Axisymmetric  contoured  nozzle  or  flexible  plate 
Figure  C-3.  Simple  relationship  for  defining  the  tunnel  wall 
boundary-layer  thickness  (u/ue  —  0.99). 
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Sym  Moo 
O  1.5 
□  2.0 


Figure  C-4.  Comparison  of  estimated  and 
measured  RE#  parameter 
in  Tunnel  A. 
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Figure  C-5.  Estimated  local  friction  coefficient  on  the  nozzle 
wall  in  the  tunnel  test  sections  over  the 
operational  Reynolds  number  ranges. 
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Sym 

O  Tunnel  A  Side  Wall 

•  Tunnel  A  Flexible  Plate 

□  Tunnel  B  and  C  Axisymmetric  Nozzle 


Mach  No. 

a.  5/6* 

Figure  C-6.  VKF  tunnel  wall  turbulent  boundary- 
layer  shape  factors. 
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b.  5/6* 

Figure  C-6.  Concluded. 
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Table  C-1  Summary  of  Tunnel  Conditions 
a.  Tunnel  A 
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Table  C-1.  Continued 
a.  Concluded. 
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Table  C-1.  Concluded 
b.  Tunnel  B 


St 

C 

illini 
■i  amber 

i 

Test  Soct.  ion 

Probe  Measurements 

Mach 

P, 

RHO, 

MU, 

■ i 

r- 

Hln, 

No. 

PT, 

TT, 

DP, 

Too  f 

psl  a 

Q, 

Lbm/ f  t3 

u 

I 

(1) 

ft-1 

PP, 

MUO 

fl-l 

psia 

OR 

or 

or 

x  10-2 

psia 

x  102 

ft/SeC 

x  I07 

x  10-H 

psia 

ML 

x  107 

X  10~5 

5 .  bS 

2G.4 

H 

15 

107.4 

1.89 

0.458 

0.0476 

2.S 

90 

0.863 

0.51 

0.801 

0.768 

6.03 

5. 

35 

0.828 

5.89 

53.7 

107 . 0 

3.8J 

0.925 

0.0960 

0.861 

1.04 

1.72 

1.61 

5.99 

5. 

35 

1.68 

b.89 

79.5 

1 

107.0 

5.64 

1.07 

0.142 

0.861 

1 . 53 

2.55 

2.35 

6.00 

N. 

D. 

2.88 

5.91 

107.1 

l] 

106.4 

7.44 

1.82 

0.189 

0.856 

2.05 

3.38 

3.30 

5.95 

5. 

37 

3.24 

5.92 

160.0 

1 

106.1 

11.00 

2.70 

0.280 

0.853 

3.05 

5.02 

4.99 

5.93 

4.80 

5.94 

212.0 

105.5 

14.28 

3.53 

0.365 

0.848 

4.00 

6.55 

6,59 

5.93 

6.36 

5.95 

273.  3 

■ 

105.2 

18.22 

4.52 

0.468 

0.846 

5.14 

8.39 

8.38 

5.95 

5.38 

8.04 

5.92 

159.0 

1,115 

-55 

139.0 

10.93 

2.68 

3.420 

7 

-Stage 

Operat 

Ion 

7.89 

115.5 

1,340 

1 

99.6 

1 .29 

0.564 

0.0351 

3,360 

o.soi 

0.525 

1.04 

0.950 

8.06 

7.23 

0.588 

7.94 

231  .6 

1 , 360 

99.9 

2.40 

1.10 

0.0673 

3,890 

0.803 

1.01 

2.03 

1.92 

8.04 

7.31 

1.08 

7.99 

462.3 

1  ,350 

98.  1 

4.77 

2.13 

0.131 

3,880 

0.788 

2.01 

3,95 

3.84 

8.04 

7.27 

2.16 

8.00 

693.0 

1 

97.8 

7.10 

3.18 

0.196 

0.786 

3.00 

5.88 

5.66 

8.07 

7.28 

3.24 

8.00 

849.4 

1 

97.8 

8.70 

3.90 

0.240 

0.7S6 

3.68 

119 

6.93 

8.07 

7.28 

3.96 

AEDC -TR-81-5 


AEDC-TR-81-5 


APPENDIX  D 

SIX-DEGREE  SHARP  NOSE  CONE  SURFACE  AND  FLOW-FIELD  PROPERTIES 

The  measurements  obtained  on  the  6-deg  sharp  nose  cone  were  normalized,  when 
appropriate,  by  the  surface  properties  given  in  this  Appendix  (see  Fig.D-1).  The  local  surface 
velocity  measurements  at  Mach  number  4  and  6  are  based  on  the  local  mean  flow-field 
measurements,  namely  the  local  wall  static  pressure,  the  pitot  pressure  just  above  the 
boundary  layer,  and  the  stream  total  temperature,  The  theoretical  estimates  of  the  local  cone 
properties  were  based  on  the  inviscid  sharp  cone  solutions  of  the  Taylor-Maccoll  equation 
(Ref.  45).  There  is  agreement  between  the  theory  and  experimental  values  of  the  local 
velocity  and  cone  surface  pressure  ratios  at  the  maximum  free-stream  Reynolds  number. 

Assuming  a  Prandtl  number  of  0.71,  a  comparison  between  the  measured  cone  surface 
temperature  and  the  estimated  value  is  reasonable  in  the  Mach  number  range  from  1 .5  to  5. 
The  lower  experimental  value  at  Mach  number  6  was  obtained  in  Tunnel  B  where  the  model 
was  not  allowed  to  reach  a  stable  (essentially  adiabatic)  condition. 

The  estimated  laminar  boundary-layer  thickness  for  the  6-deg  sharp  nose  cone  at  the 
upper  and  lower  free-stream  Reynolds  number  is  shown  in  Fig.  D-ld.  These  theoretical 
estimates  were  based  on  results  derived  from  Refs.  46  and  47.  The  corresponding 
experimental  results  at  Mach  Nos.  4  and  6  were  based  on  mean  flow  measurements  made  at 
model  station  25  (in.)  on  the  cone.  Unfortunately,  the  14  percent  disagreement  between  the 
Mach  No.  6  resulls  and  the  theory  was  probably  due  to  a  large  bias  in  the  uncertainty  in  the 
pitot  probe  displacement  measurement  used  to  define  the  boundary-layer  profile.  In  this 
study,  the  analytically  defined  parameters  were  selected  for  normalizing  cone  surface 
pressure  fluctuation  measurements. 
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Nominal  Free-Strearo  Mach  Number 

b.  Average  surface  pressure  ratio 
Figure  D-1.  6-deg  sharp  nose  cone  surface  properties. 
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•Adiabatic  Wall 

Average  Model 
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Condition  (Pr  =  0.71) 

Wall  Condition^ 

_ 1 _ 1 _ 

_ 1 _ 

c.  Average  model  wall  temperature 


Nominal  Free-Stream  Mach  Number 


d.  Boundary-layer  displacement  thickness 
Figure  D-1.  Concluded. 
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NOMENCLATURE 

A  Constant 

Aj  Test  section  cross-sectional  area,  in.2 

a  Speed  of  sound,  ft/sec 

B  Intercept  of  (SF  =  m  P  +  B) 

BW  Total  bandwidth  of  the  analysis,  Hz 

C(  j  Coefficient,  (  )/qro 

cf  Local  wall  shear  stress  coefficient 

Dj  Coefficient,  j  =  1  to  5 

DE  Effective  diameter,  DE  =  (4AT/7r)1/2,  in. 

DP  Dew  point,  °F 

e  Voltage  across  wire,  volts 

Ae(  >  Sensitivity  coefficient  of  wire  for  (  )  component,  volts 

F  Power  spectral  density  function,  units2/Hz 

Fo  Fg  Winter  and  Gaudet  (Ref.  43)  correlation  parameters 
f  Frequency,  Hz 

61*82  Coefficients  defined  by  Eqs.  (18)  and  (19),  respectively 
H  Transform,  units/Hz 

h  Real  time  signal,  units 

I  Imaginary 

k  Number  of  the  discrete  point  in  the  real  time  signal 

LI,  L2  Pitot -acoustic  probe  support  length,  in. 
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Moo 

ML 
MUO 
MU,  p 
m 

N 

NED 

Nua 

n 

nd 

P 

Pi.  P2 

Pp 

Pr 

PST 

PSW 

pi 

PT2 

P*.  PW( 


Mach  number 
Local  Mach  number 
Viscosity  based  on  TT,  lbf-sec/ft2 
Viscosity  based  on  T,  lbf-sec/ft2 
Slope  of  (SF  =  m  P  +  B) 

Number  of  discrete  points  in  a  signal  sample  (N  =  1 ,024  points) 

Normalized  energy  ratio  (see  Fig.  30) 

Nusselt  number  based  on  TT 

Number  of  the  discrete  amplitude  point  in  the  frequency  domain 
Number  of  signal  samples 

Directional  cosine  of  the  normal  to  the  sound  plane-wave  from  relative  to  the 
mean-flow  direction 

Static  pressure,  psia 

Pressure  upstream  and  downstream,  respectively,  of  a  normal  shock,  psia 
Pitot  probe  pressure,  psia 
Prandtl  number 

Stilling  chamber  pressure  fluctuation,  psi 

Side  wall  static  pressure  fluctuation,  psi 
Stilling  chamber  pressure,  psia 

Computed  toLal  pressure  downstream  of  a  normal  shock,  psia 
Wall  or  surface  pressure,  psia 
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Tunnel  window  wall  static  pressure  fluctuation,  psi 
Dynamic  pressure,  psia 
Real  number 

Resistance  of  wire  at  condition  (  ),  ohm 

Correlation  coefficient  for  components  indicated  by  subscripts 

DPS  gage  radius,  in. 

Reynolds  number  based  on  T* 

Reynolds  number  based  on  TT 
Density,  lbm/ft3 

Boundary-layer  shape  factor  r 

Scale  factor  of  DPS  gage,  psi/mv 
Time,  sec,  or  sample  length 

Temperature  (typically  subscripted  or  superscripted),  °R 
Temperature  of  wire  at  condition  (  ),  °R 
True  rms 

Stilling  chamber  or  free-stream  total  temperature,  °R 
Local  wall  temperature,  °R 
Velocity,  ft/sec 

Overheating  parameter,  Aer/Aeff 
Surface  length,  ft 

Nondimensional  rms  voltage  fluctuation,  e7Aec 
Temperature  linear  coefficient  of  resistivity,  (°C)_] 

Ratio  of  specific  heats  (7  =  1.4) 
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Af  Bandwidth  between  points  in  a  spectrum,  Hz 

A(  )  Incremental  difference  of  the  property  (  )  from  its  mean  value 

8  Boundary-layer  thickness,  ft 

6*  Boundary-layer  displacement  thickness,  ft 

e  Uncertainty  parameter 

t)  Recovery  factor 

6  Boundary-layer  momentum  thickness,  ft 

A  Integral  scale  of  parameter  (  ) 

x  Sound  mode  “amplitude,”  normalized 

a  Entropy  mode  “amplitude,”  normalized 

t  Local  shear  value,  psi,  or  vorticity  mode  “amplitude,”  normalized 

<(>  Phase  angle  or  model  roll  angle,  deg 

a>  Angular  velocity,  27rf,  radians  per  second 

SUBSCRIPTS 

2  Downstream  of  a  normal  shock 

C  Convective  quantity 

CORR  Corrected  or  input  value 

D  Discrete  or  digital  value 

DE  Effective  test  section  diameter,  2(area[es,  section/ *0 1/2 

d  Wire  diameter 

f  Pertaining  to  reference  temperature  of  32°F 

i  Local 

m  Measured 

o  First  mode 
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r  Reference  value 

ST  Stilling  chamber 

STA  Model  (or  tunnel)  station 

s  Source  (or  shock)  velocity 

t  Transition  location 

w  Pertaining  to  temperature  of  heated  wire  or  wall 

x  Based  on  surface  length 

0  Based  on  boundary-layer  momentum  thickness 

ir  Pertaining  to  sound  mode 

a  Pertaining  to  entropy  mode 

t  Pertaining  to  vorticity  mode 

w  Free  stream 

SUPERSCRIPTS 

(~)  Square  root  of  the  mean  square  value 

(')  Fluctuating  component 

*  Complex  conjugate  or  reference 

(")  Average  with  respect  to  time 
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